
'"i . av ■ - 

"A -V \i 
■■ir't-r ■ 




- s..y k r .'-:\ "•* j- 


A 
ri C£ 


r, i r v v - v ..x 

w, a <* " 


;/• 


SkAAAAy 

• - ■' ~ j 


’ r . “ , ' .' r 'i Ss, -* *-''-*'3 

.■V ,» m r*s. 1 — 

V . A -’rV> 


. < A. -rvv ' 

t ,-^ ■/'** ' 


:-. ru : ) ; * ^ :•>>;• % 

’ GODDARD BROUWER ORBIT BULLETIN 

r te 4 % ■ ;&> 

vy-A -. xa, / ■ ,y ■ 7 \ h -■■ • . - ; ^r.-- • ,. Tt '- -*£. , ; </ . ~--v A s -- < ’ 0 ;' f f^> ... • a* i , - 

, >.' ’ ->>■ >: ‘V^ a\ . *>• ■-:*'«■. --*y?r 

^ r. .-*. I" ‘* x ~ - ' - •— ;f , - -• >- • j* ■__. „ v //' -yT'V • . ■ IV^V »' ]£, “. z**'* . - '. 

yiCjs- X yr<X j -:X J $'^ ~ '■-■■ "< ■■ ^ - ' a-' -v -• yy-A&C- CtvVb'^ • >-*•' A 

s-vS: ;* =r-Vr~xx ,'^vi;,r.;. J--/A 3 A 


.r-. V<' 7 '-^ 


C; 4 ^'” : ^'r 

^ A v • ‘ y Vs 
•’«■ Jfr' v r - .O 


, - .. , 
.;V^K' ^V- 


SV ■ '■’ , > . 

■ • *' 

o ‘v- - ,- « "s -' 


'Si 

7- 


_v 


-f \„^ v -' af* .• . (:/'■ -A' •>.'. M ; , >'A', '-> v , ',;*>£ ->>Sf'M -s 

vl# ;-.^D " >>• c i -^K^' v '’’-- ' a- r ^ r 

V.S -,1 iK?oRM N A A nrS CHNICAl ; 

^ ^ „' , ' > i .; ,■■£-- ■ -s- • - -ns ' l ] 'formation service 1 

A> -y/vfV&j : 4 i’&Zi&u 4 ' x V 44 x' a - jagr^i'issr- E 



A V' 


' v"c/; , 


t L r I blU 

„.. - . GREENBELT. MARYLAND ' ' 

o v-frA.-vi-^ -•■ '-^’•' -’ 


J BULLETIN (NASA) 200 p HC | i ""- T 


N73-15869 


iuu \ 

^ J CS^L 22 C 






- '.-V?-- '■ W -• - •■■ :•;.- , 

;' r V U ;.V ,y ; '•- '" 


Unclas 

G3/30 52964_ ^ , 

-• r "' r ' k '''"^ 


?, -- A, ;> 

k - --V ^ A' ’.A 


A r V 


- — Tj-r-i — 7 ^^ — r 73 ? ! ^' "••* r 7 ' ■ •■ J K'rf 5 ^" r '-" '-' v . .'V'-^-^- 

■' - • — > -.• »- ..'- VA • ■ V /-'•C.f « ,.<••'-•< Ar k > 


• '' -■ ■'"- ;. -n A ^ - .-<: ^ it- <-■' c . ^ 
/y 'A f.-A * . y- N . - ’ y ^ ^y 


■: 

. c^ r y r y V-< .’ a 


GODDARD BROUWER ORBIT BULLETIN 


D. B. Morgan 
R. A. Gordon 


July 1971 


GODDARD SPACE FLIGHT CENTER 

Greenbelt, Maryland 
• 

I 


ACKNOWLEDGEMENT 


The authors wish to gratefully acknowledge the 
guidance and technical advice given by Mr. William 
N. Weston from the Trajectory and Dynamics Branch, 
Trajectory Analysis and Geodynamics Division. 


11 



GODDARD BROUWER ORBIT BULLETIN 


D. B. Morgan 
R. A. Gordon 

ABSTRACT 

The Goddard Brouwer Orbit Bulletin provides operational 
support for earth space research and technological missions by 
producing a tape containing pertinent spacecraft orbital infor- 
mation which is provided to a number of cities around the world 
in support of individual missions. This document presents a 
program description of the main and associated subroutines, 
and a complete description of the input, output and requirements 
of the Bulletin program. 
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GODDARD BROUWER ORBIT BULLETIN 


I. INTRODUCTION 


Goddard Brouwer Orbit Bulletin is an economical means of providing opera- 
tional support for earth space research and technological missions. The 
Bulletin routine accepts as input a set of orbital elements generated by the 
Definitive Orbit Determination System (DODS) and produces an output' tape 
containing pertinent spacecraft orbital information. This information is pro- 
vided to a number of cities around the world in support of missions such as 
International Satellite for Ionospheric Studies (ISIS). The Bulletin informa- 
tion includes the following: 

1. The mean characteristics of the orbit of the satellite at epoch 

2. Prediction space elements for use when approximate satellite posi- 
tions are needed 

3. Osculating space elements 

4. Ascending nodal crossings during a requested time period 

5. An ephemeris which furnishes the positions of the satellite at 
regular intervals 

6. Brouwer data acquisition facility parameters to be used by each 
data acquisition facility to generate its topocentric predictions for 
satellite acquisition 

Section II is a program description of the main and associated subroutines. 

A complete description of the input, output and requirements is given in 
Section III, Operating Instructions. 


1 



II. PROGRAM DESCRIPTION 


A. PURPOSE 

This program provides an economical means of disseminating pertinent space- 
craft orbital information to observing stations and other interested parties. 


B. FLOW CHARTS AND FUNCTIONAL DESCRIPTIONS 


The following pages contain flow charts and functional descriptions of the main 
routine and associated subroutines. Flow charts and corresponding descrip- 
tions are grouped alphabetically, with the main routine presented first and the 
block data last. 



ENTRY 

OR 

EXIT 



DECISION 



CONNECTOR 


OPERATION 


READ 

CARD 



SUBROUTINE 
(OR FUNCTION) 
XXXX I s the Name 
of the Subroutine 




Figure 1. Flow Chart Symbols 


OUTPUT 
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MAIN Flowchart 

































BRWORB 


Brouwer Orbit Generator 

PURPOSE 

Given a time t, referenced to some epoch, the subroutine determines a set of 
osculating elements corresponding to this time with the Brouwer Orbit Theory. 

METHOD 

Brouwer (59) made use of Von Zeipel's procedure to modify Delaunay's method 
in the development of an artificial satellite theory. The subroutine "BRWORB" 
is a faithful coding of Brouwer's formulas as they appear in Sec. 9 "Formulas 
for Computation" , Brouwer, D., "Solution of the Problem of Artificial Satellite 
Theory Without Drag," Astronomical Journal, 64 (November 1959), 378-397. , 
except for modifications made to include the perturbation (pert) tape option. 

FORMULATION 

I. Compute Brouwer epoch elements corrections: 

1. Without Pert tape 



2. With Pert tape 

a. For £" calculation 

a" - a Q + Aa 


6 



e" - e Q + Ae 
i" = i Q + Ai 



b. For g" and h" calculation 



c. Correction to 2 0 , g 0 , and h 0 (epoch angular elements) 


e 0 di 


B 0 

h„ 


e 0 + a£ 


- g 0 + A g 


h Q + Ah 


2" 


at (t = t Q ) 


II. Calculation of abbreviated notation to simplify formulas: 


V = 


7 2 


A 

k. 


>■'2 


9 = cos i" 


. " 2 


7 4 


a 


"4 


72 


72 


74 


74 

•v.8 
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in. Compute the first time derivative of the secular terms: 

1. Mean Mean anomaly derivative, Anomalistic Mean Motion and period; 


e = 


di" 

dt “ n o t 


n oV 7 2 


-| (39 2 - 1) + 7 2 ' [25t7 2 + 16 7] - 15 


+ (30 - 96t, ~ 90 r) 2 ) 6 2 + (105 + 144t, + 25t, 2 ) 0 4 ] 


+ 7 4 ' e 2 (3 3O0 2 + 35 4 ) 1 


J 


• 277 

1 o D = n o + 1 > n = kl 0 D, P = — , i 0 = n 0 ‘ + i at (t = t Q ) 


2. Mean Argument of Perigee derivative; 


• _ dg 

g - 


& = no {y 


( 5(9 2 - 1) + -^2 y 2 ' f 25 77 2 + 247, - 35 


+ (90 - 192 7 ] - 126 t, 2 ) 0 2 + (385 + 360t, + 45 t, 2 ) <9 4 ] 
+ 7 4 ' [21 - 9 T7 2 + (126t, 2 - 270 ) 6 2 + (385 - 189t, 2 ) 9 4 ] \ 


3. Mean longitude of ascending node derivative; 


dh M r 

3 

- 

dt _ n 0 j?2 

.8 7 

2 [(9t , 2 + 1277 - 5) 9 + (-35 - 36t, - 5t, 2 ) 0 3 ] - 39 


+ -f-7 4 ' ( 5 - V) 6 ( 3 - 79 2 ) 
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IV. Compute constants for long-period terms; 


3 


(1 - 50 2 * )- 


Compute CP 1 - CP 15 

CP; = 406> 4 (1 - 502)' 


i 8 7 2 


y' v 2 (1 - H£ 2 - £P 2 ' ) 


e p 2 - ~~ v 2 [1 - 302 - 80 4 (1 - 50 2 )- 1 ] 

^2 


1 ^3 

ep 3 - j— 7,2 Sin i" 


eP 4 ' = [1 - 90 2 - 240 4 (1 - 50 2 )- 1 ] 


5 y 

CP 4 s 64 772 CP 4 sin in 

7 2 


CP' = 1 - 502 - 160 4 (1 - 50 2 ) - 1 


£P= = -^T — e" 2 -n 2 sin i" £P e 


s 384 


gP 6 ' - CPj (1 - 50 2 ) - 
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GP g = 3 + 16 e 2 (1 - 56> 2 ) -1 + £P 6 ' 

GP ? = 4 + 3e" 2 


35 7s 

GP g = e" 3 0 sin i" ( 5 + 320 2 (1 - 5# 2 )" 

^2 


-1 


GP„ 


s" e 


sin l 


35 y $ 

^ p io = 1152 7 e * P 9 £P s 

y 2 


CP 


sin l 


11 


CP 12 a 2 + e" 2 


CP U ■ (2 + 3e" 2 ) 0 2 


CP 14 - 8 (2 + 5e" 2 ) 0 4 (1 - 50 2 )' 


-l 


CP 


15 


T] 2 tan i" 


Compute Aj - Ajj 


Aj = e" (CP l - CP 2 ) 

A 2 = CP 3 + (4 + 3e" 2 ) CP 4 


A 3 = n (CPj - CP 2 ) 


A 4 = -J- [Gp s + (4 + 9 e" 2 ) CP 4 ] 


+ 2CP 6 ') 
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A 5 = ir lp ! 

A, = i y‘ («P ia - nSP 13 - 5SP 14 - 10e" 2 0 2 8P e ' ) 

5 74 

+ 24 ~T ( £P 12 - 3CP 13 - e p !4 - 2e" 2 2Pg ) 
72 

a 7 = i — (* p n - 9 fip 9 > + 4 — [C 7 ? 2 «p u - ^ *p 9 > *:p 7 

y 2 • r 2 


+ 


sin i" (26 + 9e" 2 )] CP 4 ' - 


_15 73 

32 ri 


e" e 2 sin i ep ? gp 6 


35 7s_ 

7 2 


"8 ^ U52 — [e" sin i" (3 + 2e" 2 ) - e" 2 0CP 9 ] £P S ' + 0»P 8 


1 5 74 

A g = - 7 2 e" 2 6 [11 + 8O0 2 (1 - 50 2 )" 1 + 5 C P 6 ‘ ]' + e" 2 0E'P 6 

8 lz . 


A 10 = £P 9 


1 73 5 7s 

4 , + 

7 2 7 2 


j?p Pp ' 

64 , X1 7 Xi 4 


A 11 «PlO + CP 8 


V. Compute constants for short-period terms included: 


Compute SP! - SP 6 


SP = — — 
1 2e" 

SP 3 - "2 7 2 ' 

sp 2 s r 2 ( 1 - ^ 2 ) 

SP 4 = 2 (- 1 + 30 2 ) 
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sp 5 ■ 6 (- 1 + so 2 ) sp 6 = esY 3 {\-d 2 ) A 

VI. Call DRAG Subroutine to compute Al drag at Observation time t. 

m 3 

= E E N p (t - t q ) p 

q=0 P=2 

where 

m = 0, 1, 2, . . . , 19 

VII. Compute Secular Terms: 

1. £" = mean mean anomaly; 

fiodgL = ( e o + AC drag) £ o d g pt = mod (£ 0 d£ + A£ drag , 2-77) 

£" = mod (n 0 't, 277 ) + mod(£t, 2tt) + mod(£ 0 d£ + AC drag> 2tt) 

2. g" = mean argument of Perigee; 

g" = mod (gt + g Q , 2 t7) 

3. h" = mean longitude of Ascending node; 


h" - mod(ht + h 0 , 2n) 

VIII. Test for Critical Inclination: 

Ai = |i" - i c | 

where 


if 


i c = 63.43° 

Ai < 1.5 
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then 


S 1 e = S 1 I = l'=g'-h' 

IX. Compute Long-Period terms: 

1. £' = mean anomaly; 

£' = £" + A 3 sin 2g" - A 4 cos g" + A g cos 3g" , mod (£', 2n r) 

»! = «' - «„t + e 0 + A£ drag , mod (£j, 2 tt) 

2. g' = Argument of Perigee; 

g' = g" + A 6 sin 2g" + A ? cos g" + Ag cos 3g" , mod (g' t 2rr) 

3. h' = longitude of Ascending node; 
h' = h " + Ag sin 2g" + A 10 cos g" - Aj 3 cos 3g" , mod (h 1 , 2rr ) 

4. Call KEPLR1 Subroutine to determine E' and compute f ' from 

f ' = tan -1 

it 

r ' ” 1 - e" cos E' 

X. Compute Short -Period terms included: 

Compute Bj - B 6 

Bj - 7 2 (-1 + 30 2 ) - t?- 3 ) + 3(1 - 0 2 ) cos ( 2 g ' + 2 f ’) 

B 2 = 3e" cos ( 2g ' + f ') + e" cos (2g‘ + 3f ') 


/l - e" 2 sin (E'> 
cos E' - e" 
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+ 


a 


B, 


B. 



SP 4 (B 3 


r' 

+ 1) sin f ' 


+ 3(i - e 2 ) 


(-B 3 + 1) sin ( 2 g ‘ + f ') + 


( B 3 + I") 


sin (2g' + 3 f ’) 


B s = f ' - 1' + e" sin f ' 

B 6 = 3 sin ( 2 g ' + 2 f ’) + 3e" sin(2g' + 2 f ' ) 

+ 3e" sin (2g‘ + f ') + e" sin (2g‘ + 3f ’) 


XI. Compute Osculating Elements: 

1. Compute a (semi-major axis) 

a = a" (1 + B x ) 

2. Compute e (eccentricity) 

e = e" + S x e + SP X (b x - 7 2 tT 4 cos (2g' + 2f) 3(1 - 02) - SP 2 B 2 ) 

3. Compute i (inclination) 

i = i" + Sji + SP 6 [3 cos (2g‘ + 2f ') + B 2 ] , mod (i, 2 tt) 

4. Compute g (argument of Perigee) 

BP 3 

g - g'+SP 1 SP 3 B 4 + — [SP 5 B s + (3 - 50 2 ) B 6 ] , mod(g, 2 tt) 

5. Compute h (longitude of ascending node) 

h - h - 0 SP 3 (6B S - B 6 ) , mod (h, 2tt) 
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6. Compute 1 (mean anomaly) 

8 = 8' - 7 ? SP 1 SP 3 B 4 , mod (8, 2tt) 

7. Call KEPLR1, to compute E (eccentric anomaly) 

8. Compute f (true anomaly) 


f 


tan 


-l 


yi - e 2 Sin E 
cos E - e 


XII. Compute Position and Velocity Vectors: 

Call the UVUK Routine to perform the mapping, (osculating keplerian 
elements to rectangular cartesian) 

a, e, i, g, h, 8 - x, y, z, x, y, z 


CALLING SEQUENCE 

CALL BRWORB (TO, T, DRG, NQ, F, EA, R, DR, RMAG, DRMAG, N, PD, 
PASS, K, SATID) 


INPUT/OUTPUT 


Arguments 


I/O 

Variable 

Description 

I 

TO 

Epoch time and date 

I 

T 

Observation time 

I 

DRG(60) 

DRG(l) 1 

j, r t 0 , 1 1, . . . 1 19 

DRG(20) J 
DRG(21) 'j 

I f N 2 , o» N 2> i» . . . N 2 lg 

DRG(40) J 
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Arguments (continued) 


I/O 

Variable 

Description 



DRG(41) 'l 

1 f ^3,0’ ^3, 1* * * * *» ^3, 19 

DRG(60) J 

I 

NQ 

Number of drag inputs 

o 

F 

True Anomaly 

o 

EA 

Eccentric Anomaly 

o 

R(3) 

x, y, z Satellite Position Vector 

o 

DR (3) 

x, y, z Satellite Velocity Vector 

o 

RMAG 

r — magnitude of Position Vector 

o 

DRMAG 

v — magnitude of Velocity Vector 

o 

N 

Anomalistic mean motion 

o 

PD 

Anomalistic Period 

I 

PASS 

PASS = 1 Compute constants (at t^ needed in computa- 
tion of osculating elements 
PASS = 2 Update osculating elements to observation 
time t 

I 

K - ^GM 

Gravitational Constant (length) 3 7 2 /time 

I 

SATID(ll) 

SATID(l) = Satellite identification number 

SATED (2) = reference year 

SATID(ll) = day count of reference date 


Common 


I/O 

Block 

Variable 

I 

BPOOL 

TABLE(12) - TABLE(15) = J 2 , J 3 , J 4 , J 5 
TABLE (16) = Deg. /rad. 

TABLE (22) = r e 
TABLE (31) = TOL 
TABLE (41) = n 

TABLE (61) - TABLE (64) = K 2 , K 3 , K 4 , K 5 

O 

O 

SECPRM 

DOTELE 

DPELE(6) = a", e", i", g", h", 1" 

LODOT =J 0 

LDOT = 8 

GDOT = g 

HDOT = h 
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Common 


I/O 

Block 

Variable 

o 

LPPRM 

DELIE = Sje 
DEL1I = 5ji 
LI = fi, 

LP = 8' 

- 


GP = g' 
HP = h' 

o 

OSCELE 

ORBPRM(6) = a, e, i, g, h, £ 

o 

ETAP 

V 3 , V 6 , V 4 

o 

THETA 

M1P3T2 = 30 2 - 1, THETA = 0 = cos i 

o 

GMPR 

72 

o 

UVPQ 

U, V, P, Q 

o 

DGPRM 

LODGL = + AE(j ra g 

o 

DAFPRM 

LDAF(l) = Aj 
LDAF(2) = A 2 
LDAF(3) = -EP 5 
LDAF(4) = A 3 
LDAF(5) = -A 4 
LDAF(6) = A 5 
LDAF(7) = A 6 
LDAF(8) = A 7 
LDAF(9) = A 8 
LDAF(IO) = A 9 
LDAF(ll) = A 10 
LDAF(12) = -A n 
LDAF(13) = £P 15 

I 

PIND 

NPT 

o 

PERTL 

LODGPT = £ 0 + A£ + AE^g 

E 0 

o 

DELKEP 

DKEP(6) = Aa, Ae, Ai, Ag, Ah, A£ 

o 

PRTKEP 

PKEP(3) = g 0 + Ag, h 0 + Ah, £ 0 + A£ 

o 

NOD 

LOD = £ 0 D 


1 — 1 

CALLED BY 

CALLS 

DAF 

PQ 

ELEMLD 

PE RTF 0 

PREDS 

DRAG 

SPACE L 

KEPLR1 

GTRACE 

REDUCE 

NODALX 

UVUK 

NSPT 
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BRWORB Flowchart 






























































CHANPL 


Change of Constants 


PURPOSE 

To change any constant that is in the BLOCK DATA or POOL Subroutine. 


METHOD 

Constants from the BLOCK DATA or POOL subroutine are changed according 
to the values on the change of constants cards. Each card permits change of 
one to three constants. The first forty constants, which are in the BLOCK DATA, 
can be changed by the first change of constants card(s). The remaining forty 
constants, which are defined in the POOL subroutine, can be changed by the 
second change of constants card(s). All eighty constants are stored in COMMON 
BPOOL. 


CALLING SEQUENCE 
CALL CHANPL 


INPUT/OUTPUT 


Common 


I/O 

Block 

Variable 

Description 

I/O 

BPOOL 

TABLE (80) 

Frequently used constants that are 
previously set in the program but may 
be changed by the change of constants 
cards. 


CALLED BY 
MAIN 
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CKSUM 


Check Sum 


PURPOSE 

To sum the digits of a line of data to modulo ten. 


METHOD 

A data line in SPACEL contains ten numbers which will be summed modulo 
ten in CKSUM. 


CALLING SEQUENCE 


CALL CKSUM (NCK, NX) 

INPUT/OUTPUT 

Arguments 


I/O 

Variable 

Description 

I 

0 

NCK(IO) 

NX 

Array of ten numbers from a line of data 
Sum of all digits from NCK modulo ten 


CALLED BY 
SPACEL 
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DAF 


Data Acquisition Facility Parameters 


PURPOSE 

To compute Brouwer parameters to be used by each data acquisition facility to 
generate its topocenfrric predictions for satellite acquisition. 


METHOD 

The following quantities are computed by Brouwer Satellite Theory Orbit Gen- 
erator (BRWORB). They are defined here using Brouwer's Terminology. 


LDOT 


1 + ~2 y' 2 V (“ 1 + 30 2 ) + y' 2 2 17 [- 15 + I 677 + 2Srj 2 


+ (30 - 96 7) - 90t7 2 ) 6 2 + (105 + 144t] + 25t? 2 ) 0 4 ] 


+ -jg- y\ rj e " 2 [3 - 30 <9 2 + 35 0 4 ] 


GDOT 




+ (90 - 1927] - 126t7 2 ) 6 2 + (385 + 360t} + 45t] 2 ) 6 4 ] 


+ y\ [21 - 9t ? 2 + (-270 + 126 77 2 ) 6 2 + (385 - 189 t? 2 ) 0 4 ] 
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HDOT = 


LDAFi = 


LDAF 2 = 


LDAF 3 = 


ldaf 4 = 


- 3 y' 2 d + y' 2 [(“ 5 + 12t] + 9 t 7 2 ) 0 + (- 35 - 36r} - Srj 2 ) 0 3 ] 

+ 4 7; (5- St? 2 )^ (3- 70*) - 


72 e" T7 2 [1 - 1102 _ 4004 (1 _ 502)-!] 


- — e " 772 [ i - 302 - 8 0 4 ( i - 50 2 ) -1 ] 


7 2 



sin I 


+ 


_ 5 ^ 

64 


7; 


7 2 V 


sin I" (4 + 3 e" 2 ) [1 - 9 0 2 - 24 0 4 (1 - 50 2 )- 1 ] 


- “ e" 2 77 2 sin I" [1 - 502 - 16 0* (1 - 50 2 )- 1 ] 

y 2 


8 7 2 


7 ' 77 3 [1 - 1102 _ 4 Q 0 4 (1 - 50 2 )- 1 ] 


“ 12 ~y 773 Cl ‘ 302 ” 804 C 1 " 50 2 ) -1 ] j 
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ldaf 5 = 


ldaf 6 = 


ldaf 7 = 


LDAFg = 


4 7; e" 


sin I" 


5 

64 


3 

sin I" (4 + 9e" 2 ) [l - 9 9 2 - 24 0 4 (1 - 50 2 )' 1 ] • 

72 e 


35 y's 

384 v 
7 2 


?7 3 e' 


sin I 


[1 - 50 2 - 16 0 4 (1 - 502)-i] 



7j [+ (2 + e" 2 ) - 11 (2 + 3e" 2 ) 6 2 - 40 ( 2 + 5e" 2 ) 6> 4 ( 1 - 56> 2 )-i 


- 400 e" 2 £? 6 (1 


50 2 )- 2 ] 


_5_ y\ 
24 7i 


[2 + e" 2 


3(2+ 3e" 2 ) 6 2 


-8(1+ 5e" 2 ) 0 4 (1 - 5 0 2 )' 1 - 80 e" 2 0 6 (1 


5# 2 r 2 i 


1 

73 

( sin I" 

e" 02 \ 

l 4 

y'i 

( e" 

sin r) 


_5_ 2i 
64 7i 



e" g 2 \ 
sin I"/ 


(4+ 3e" 2 ) + e" sin I" (26 + 9e" 2 ) 


[1 - 9 e 2 


- 240 4 (1 - 50 2 )” 1 ] 


32 y 1 
7 2 


e" 0 2 sin I" (4 + 3e" 2 ) 


[3 + 166 2 (1 - 50 2 )- 1 + 40 6> 4 (1 - 50 2 )- 2 ] 
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LDAF 9 


r_ 35 

t" 1152 y' 2 


e" sin I" (3 + 2 e" 2 ) 


e " 3 0 2 ~ 
sin I" 


[1 - 50 2 


- 160 4 (1 ~ 5 9 2 )- 1 ] + 


35 y's 

576 y 

y 2 


e " 3 0 2 sin I" [5 + 3202 (1 


502)-l 


+ 8O0 4 (1 - 56>2)-2]j 


LDAF 10 = \~ Y ^2 e " 2 0 £ U + 8002 ( 1 " 502 ) _1 + 20004 t 1 ‘ 502 )" 2] 


+ 


5 

12 


Zl 

72 


e " 2 e 


[3 + 1602 (i - 50 2 )" 1 + 4O0 4 (1 - 50 2 )-2]| 


LDAF n 


1 Zl e " 0 

4 y' 2 sin I 


_ 5 _ y± e" e 

54 y 2 sin I" 


(4 + 3e" 2 ) [1 - 902 


- 240 4 (1 - 50 2 )- 1 ] 


15 7s_ 

32 -y' 
y 2 


e" 0 sin I" (4 + 3 e" 2 ) 


[3 + 160 2 (1 _ 502^1 


+ 4O0 4 (1 - 50 2)" 
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LDAF 


35 ^5 e" 3 6 


[1 - 56 2 - 160 4 (1 - 5# 2 )- 1 ] 


12 


1152 


7i 


sin I" 


35 y * 

576 


— 7 e" 3 6 sin I" [5 + 320 2 (1 - 5 6 2 )- 1 + 8O0 4 (1- 50 2 )“ 2 ] 


7 , 


LDAF 


13 


rj 2 tan I" 


where 

e" = eccentricity 
I " = inclination 
6 = cos i 0 

i 0 = inclination at epoch 
V = ]/ 1 ~ e " 2 

y' n = yjv 2n 

y \ n - K n /a n 

a = semimajor axis 

K = Brouwer's representation of the harmonics of the earth's potential 

n 


CALLING SEQUENCE 

CALL DAF (TO, DRG, JDTO, ETIME, SATID, ELEMO, NQ) 
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INPUT/OUTPUT 


Arguments 


I/O 

Variable 


Description 

I 

TO 

Epoch date and time in CUT 

I 

DRG(60) 

DRG(l) 1 
1 

DRG(20) J 
DRG(21) 1 

i 

DRG(40) J 
DRG(41) ' 

1 

DRG(60) 

' *0* 1 1* • . . 1 1 9 
► N 2j o» N 2( i» . . • N 2> i9 

1" N 3, 0» N 3. 1 N 3, 19 

I 

JDTO 

No. of days from reference to epoch 

I 

ETIME 

Epoch hours, minutes, and seconds converted to 
seconds 

I 

SATID(ll) 

SATID(l) 
SAT ID (2) 
SATID(ll) 

= satellite identification number 
= reference year 
= day count of reference date 

I 

ELEM0(6) 

Orbital elements (a, e, i, u, fl, M) 

I 

NQ 

Number of drag inputs 


Common 


I/O 

Block 

Variable 

I 

BPOOL 

TABLE (2) = KMCUL 
TABLE (24) = BK 
TABLE(35) = MINDAY 
TABLE (41) = MU 
TABLE (49) = w e 
TABLE(59) = MINCUT 
TABLE(61) - K2 

I 

DAFPRM 

LDAF(13) 

I 

DOTELE 

LDOT 

GDOT 

HDOT 

I 

THETAP 

THETA 

I 

RADIAN 

AMBDA - 
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TAPE OUTPUT 


The output from the Brouwer DAF parameters function is written on tape. Two 
lines of identifying information precede the six lines of output parameters. 

Line 1 — Title 

Col. 1 Blank 

Col. 2-24 'Brouwer DAF Parameters' 

Col. 25-80 Blank 

Line 2 — DAF Parameters Bequest Card Printout 


Col. No. 

Format 



Description 

1 


Blank 



2-8 

XXXXX 

Satellite number 

9 


Blank 



10-11 

XX 

Year 

> 


12 


Blank 



13-14 

XX 

Month 


► Date predictions begin 

15 


Blank 



16-17 

XX 

Day 

J 


18 


Blank 



19-20 

XX 

Hour 

*> 


21 


Blank 


• Time predictions begin 

22-23 

XX 

Minutes 

J 



If the satellite number of the Brouwer DAF request card agrees with the input 
satellite number then Col. 24-80 are blank. If not — Col. 30-52 of line >2 con- 
tain **ERROR IN DAF SAT ID**. 

Lines 3-8 contain the requested output parameters in the following floating 
point decimal form: 


Col. No. 

Format 

1-15 

SXXXXXXXXXXDSXX 

16 

Blank 

17-31 

SXXXXXXXXXXDSXX 

32 

Blank 

33-47 

SXXXXXXXXXXDSXX 

48 

Blank 
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Col. No. 

Format 

49-63 

SXXXXXXXXXXDSXX 

64 

Blank 

65-79 

SXXXXXXXXXXDSXX 

80 

Blank 

81-95 

SXXXXXXXXXXDSXX 

96 

Blank 

97-111 

SXXXXXXXXXXDSXX 


The order Of the parameters is as follows: 

Line 3 - LDAF 4 /2 u 
LDAF s /2 77 
LDAF /2 77 

O 

ldaf 7 

LDAFg 

LDAF 9 

LDAF 10 

Line 4 — LDAF t x 
LDAF 4 2 
LDAF t 

ldaf 2 

ldaf 3 

cos I" 
e" 

Line 5 - GDOT/13.4472 

HDOT/13.4472 
So 

K 

i o'/2” 

I "/2 77 
\( t pi )/2 77 
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Line 6 — - LDAF 13 /27r 

d 
h 
m 
a 

N/277 ■ (13.4472) 2 
n 0 + Sj/ 2 77 (13.4472) 


Line 7 — K 2 

K 3 

" K 8 

1/ 4 7 2 ' 

K g /2 W 

K, 

K 4 

Line 8 — -K g 

-k 6 

7]/4e" 7 2 ' 

t . - t. 

pi 0 

CALLED BY 
MAIN 


CALLS 

DREFOD 

JDSCUT 

BRWORB 

REDUCE 

DAYCT 
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DAF Flowchart 



































DATANO 

Double Precision Arctangent (Y/X) 


PURPOSE 

To compute a value for the arctangent between 0 and 2 PI where the tangent is 
defined by the two input arguments as ARG1/ARG2. 


CALLING SEQUENCE 

DATANO (ARG1, ARG2) 
Note that DATANO is a function. 


CALLED BY 
PRINT 



DATANO Flowchart 
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DATE 


Calendar Date 


PURPOSE 

To convert year and day count (number of days from January 0 of given year to 
a given date) to year* month and day. 


CALLING SEQUENCE 

CALL DATE (IYR, IDY, IY, IM, ID) 

INPUT/OUTPUT 


Arguments 


I/O 

Variable 

Description 

I 

IYR 

Year 

I 

IDY 

Number of days from January 0 of year IYR 

o 

IY 

Year 

0 

IM 

Month 

o 

ID 

Day 


CALLED BY 

MAIN 

SATOR 

TIMETB 
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DATE Flowchart (continued) 









DAYCT 


Day Count 


PURPOSE 

To convert year, month and day to the number of days from January 0 of a given 
year to the given date. 


CALLING SEQUENCE 

CALL DAYCT (IY, IM, ID, IDAYS) 

INPUT/OUTPUT 

Arguments 


I/O 

Variable 

Description 

I 

IY 

Year 

I 

IM 

Month 

I 

ID 

Day 

o 

IDAYS 

Number of days from January 0 of given year to 



given date 


CALLED BY 

MAIN 

DAF 

DREFOD 
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DAYCT Flowchart 














DMSTOR 


Degrees, Minutes, Seconds to Radians 


PURPOSE 

To convert degrees, minutes and seconds to radians. 


METHOD 

Radians = ((sec/60 + minutes)/60 + deg)/degrees per rad. 


CALLING SEQUENCE 

CALL DMSTOR (DEG, FM, SS, RAD) 


INPUT/OUTPUT 

Arguments 


I/O 

Variable 

Description 

I 

DEG 

Degrees 

I 

FM 

Minutes 

I 

SS 

Seconds 

o 

RAD 

Radians 


Common 


I/O 

Block 

Variable 

I 

BPOOL 

TABLE (16) = DEGRAD (deg/ rad) 


CALLED BY 

MAIN 

POOL 
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DRAG 


Compute Drag 


PURPOSE 

To compute A£ drag which provide corrections for the Brouwer Orbit Generator. 
METHOD 

m 3 

A£ drag = L L N P, q (t - t q ) P 
q=0 p=2 


where 

m = 0, 1, 2, 3, . . 19 

t = observation time 1 
* 

t q = drag time 

N p q = drag parameters 

CALLING SEQUENCE 

CALL DRAG (DRG, PI2, DRAGL, TO, T, KMULT, NO) 


INPUT/OUTPUT 

Arguments 


I/O 

Variable 

Description 

I 

DRG(60) 

DRG(l) 1 

| r t 0> tj, 1 19 

DRG(20) J 
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Arguments (continued) 


Vo 

Variable 

Description 

i 

DRG(60) 

DRG(21) ' 

^2, 0’ N 2, 1’ ’ " ° ^2, 19 



DRG(40) J 
DRG(41) ' 

1 

" ^3 i O’ ^3, 1’ • • •» ^3, 19 



DRG(60) J 


i 

PI2 

2-n radians 

o 

DRAGL 

AL 


I 

TO 

Epoch time 

I 

T 

Observation time - TO 

I 

KMULT 

K multiplier 

I 

NQ 

Number of drag inputs 


CALLED BY 
BRWORB 

CALLS 

REDUCE 
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DRAGLD 


Drag Load 


PURPOSE 

To load N q) data for the computation of AL drag. 


CALLING SEQUENCE 

CALL DRAGLD (SAT3D, DRAGDT, DRG, NQ, NERROR) 


INPUT/OUTPUT 


Arguments 


I/O 


O 


O 


O 


Variable 


SATID(ll) 


DRAGDT (40) 


DRG(60) 


NQ 


Description 


SATID(l) = reference satellite ID number 
SATID(2) = reference year 
SATID(ll) = day count of reference date 
DRAGDT(l) 


DRAGDT (3) 

i 

DRAGDT (30) 
DRAGDT(2) 

DRAGDT (4) 

1 

DRAGDT(40) 

DRG(l) 

1 

DRG(20) 

DRG(21) 

DRG(40) 

DRG(41) 

i 

DRG(60) 

Number of drag inputs 


packed drag date 

packed drag time 

1 0 > 1 1 , . . 1 19 

^2, 0 ’ ^2 , 1 ’ * * •’ ^ 2, 19 

^3, 0> ^3, 1 » • * •» ^3, 19 
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Arguments (continued) 


I/O 

Variable 

Description 

o 

NERROR 

Error indicator 
= 0 no error 

= -1 wrong Sat. ID number on drag card 


CALLED BY 
MAIN 

CALLS 

DREFOD 

JDSCUT 
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DREFOD 


Day Count from Reference Date to Observation 


PURPOSE 

To compute the number of days from the reference date to the observation date. 


CALLING SEQUENCE 

CALL DREFOD (IYREF, IDC, IOY, IOM, IOD, IDAYS) 


INPUT/OUTPUT 

Arguments 


I/O 

Variable 

Description 

I 

IYREF 

Year of reference date 

I 

IDC 

Number of days from January 0 of the year to the 
day of reference 

I 

IOY 

Year of observation date 

I 

IOM 

Month of observation date 

I 

IOD 

Day of observation date 

o 

IDAYS 

Number of days from the reference date to the 
observation date 


CALLED BY 

MAIN 

DAF 

DRAGLD 

ELEMLD 

PERTFO 

PRINT 

TIMETB 

WMAPLD 

CALLS 

DAYCT 
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ELCONO 

Elements Conversion 


PURPOSE 

• To convert position and velocity vectors to Keplerian elements. 

• To convert Keplerian elements to position and velocity vectors. 


METHOD 

ELOSC — Convert position and velocity vectors to osculating elements. 

Let r be the magnitude of the position vector r = (x, y, z) and v be the magnitude 
of the velocity vector r = (x, y, z). 

Compute 

h. = (h x , h y , h 2 ) = r x L . 

Compute the semi-major axis of the orbit: 

Mr 

a - . 

2 / j . - rv 2 

Compute the inclination angle: 


where 


tan 


-l 


(h; 




h *) 
y ' 


0 < i <7 t . 
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Compute the eccentricity: 


e 


fj.a 


- K 2 


fj. a 


where 


h 2 - h 2 +h 2 + h 2 . 

x y z 


Compute longitude of ascending node: 


Q 


tan 1 



for i / 0 . 


where 


0 < Q < 2rr 


and 


fl = 0 for i - 0 . 


Compute argument of perigee: 

co - u - v , for e / 0 , 


where: 


V 


tan -1 


~ h(r • f)~ 

. h 2 - fir . 


u 


tan 1 


z 

sin i (x cos Q + y sin Q) 


for i / 0 , 
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u 


for i 


0 , 



and 


0 _< CO < 2 TT . 

co - 0 , for e = 0 . 


Compute mean anomaly: 

M = E - e sin E, for e / 0, 


where 


E 


2 tan 1 


/ 1 - e V* sin v 
\l + e/ (1 + cos v)_ 


and 


0 < M < 277 . 

M = u , for e = 0 . 


ELIRV — Convert osculating elements to position and velocity vectors. 

Call subroutine KEPLR1 to compute the eccentric anomaly, E, given the mean 
anomaly, M, and the eccentricity, e. 

Compute the true anomaly, v : 


v 


2 tan 1 


/ 1 + e \ ^ sin E 
\ 1 - e / (1 + cos E) 


Compute the position magnitude, r: 

r = a ( 1 - e cos E) . 
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Compute the radial and horizontal components of r, V r and V p , respectively: 


V 


r 



(e sin E) 



e 


2 


Compute the position and velocity vectors, r = (x, y, z) and r = (x, y, z), 
respectively: 


x = r [cos Cl cos (a) + v) - sin Q cos i sin (co + v)] 
y = r [sin Cl cos (co + v) + cos fi cos i sin (a> + v)] 
z = r sin i sin (co + v ) 


x 


y 


Vr 

r 



r 



r 


x - V p [cos Cl sin (a > + v) + sin Cl cos i cos (co + v)] 


y + V p [-sin Q sin (co + v) + cos fl cos i cos (co + v)] 


z + V p sin i cos (co + v) . 


CALLING SEQUENCE 

Subroutine ELCONO is accessed through one of its entry points, ELOSC or ELIRV. 

CALL ELOSC (INPUT, OUTPUT); convert position and velocity vectors 
(x, y, z, x, y, z) to osculating elements (a, e, i, co, Cl, M). 

CALL ELIRV (INPUT, OUTPUT, IERR); convert osculating elements 
(a, e, i, co. Cl, M) to position and velocity vectors (x, y, z, x, y, z). 
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INPUT/OUTPUT 


Arguments 


I/O 

Variable 

Description 

From entry point ELOSC 

I 

INPUT(6) 

Position and velocity vectors (x, y, z, x, y, z) 

O 

OUTPUT (6) 

Osculating elements (a, e, i, a>, Q, M) 

From entry point ELIRV 

I 

INPUT (6) 

Osculating elements 

O 

OUTPUT (6) 

Position and velocity vectors 

O 

IERR 

Error return from KEPLR1 
= 0, convergence of eccentric anomaly 
> 0, no convergence 


Common 


I/O 

Block 

Variable 

I 

BPOOL 

TABLE (34) = PI 
TABLE (4) = MU 


CALLED BY 
ELEMLD 


CALLS 

VECTOR 

VMAG 

VCROSS 

VDOT 

KEPLR1 
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ELOSC Flowchart 
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ELOSC Flowchart (continued) 
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ELOSC Flowchart (continued) 












ELEMLD 
Elements Load 


PURPOSE 

To load the epoch and element data. 


CALLING SEQUENCE 

CALL ELEMLD (SATID, EPOCH, TO, ETIME, JDTO, ELEMO, JDTO, 
ELEMO, OSCO, PV, OUTPUT, NERROR, XINPUT, 

IF LAG, DRG, NQ) 


INPUT/OUTPUT 


Arguments 


I/O 

Variable 

Description 

I 

SATID (11) 

SATID(l) = reference satellite identification number 

SATDD(2) = reference year 

SATID(ll) = day count of reference date 

o 

EPOCH(IO) 

EPOCH(l) = epoch satellite identification number 
EPOCH(2) = year of epoch 
EPOCH(3) = month of epoch 
EPOCH(4) = day of epoch 
EPOCH{5) = hours of epoch 
EPOCH (6) = minutes of epoch 
EPOCH(7) = seconds of epoch 
EPOCH(8) = type of input epoch elements code 
= 1, position and velocity vectors 
= 2, osculating elements 
EPOCH(9) = pass number 
EPOCH(IO) = perturbation indicator 
= 0, no perturbation 
= 1, use perturbation tape 

o 

TO 

Epoch date and time in Canonical Unit of Time 

o 

ETIME 

Epoch time (hr, min, sec) in seconds 
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Arguments (continued) 


I/O 

Variable 

Description 

o 

JDTO 

Number of days from the date of reference to the 
date of epoch 

o 

ELEM0(6) 

Epoch elements in CUL 

o 

XINPUT(6) 

Epoch elements in KM 

(ELEM0(6) and XINPUT(6) may be Keplerian 
elements — a, e, i, w, n, M, or position and 
velocity vectors, x, y, z, x, y, z, according to the 
elements input type) 

o 

OSCO(6) 

Brouwer osculating elements at epoch 

0 

PV(6) 

Position and velocity vectors at epoch 

o 

OUTPUT (6) 

Converted ELEM0 elements 

0 

NERROR 

Error indicator 
= 0 no error 

> 0 elements sat. ID does not agree with reference 
sat. ED 

o 

IFLAG 

Elements unit indicator 
= 0 input elements in CUL 
= 1 input elements in Km 

I 

DRG(60) 

Drag parameters table 

I 

NQ 

Number of drag inputs 


Common 


I/O 

Block 

Variable 

0 

SECPRM 

DPELE(6) = Brouwer input mean elements 

I 

OSCELE 

ORBPRM(6) = Brouwer osculating elements 

I 

BPOOL 

TABLE (24) = BK 



TABLE(34) = PI 



TABLE (41) = MU 

o 

PERTL 

EAO = eccentric anomaly at epoch 

o 

PIND 

NPT = pertape logical number 


CALLED BY 
MAIN 
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CALLS 


DREFOD 

JDSCUT 

ELCONO 

ELIRV 

ELOSC 

BRWORB 
























GTRACE 


PURPOSE 

Routine to generate a one orbit ephemeris by advancing the satellite with equal 
intervals in geodetic latitude. 


METHOD 

Given a value for the geodetic latitude we wish to determine the corresponding 
time. A two body analysis leads to the following computational scheme. 

Given: <p s (geodetic latitude) 

Compute geocentric latitude from 

<p g - tan -1 [(1 ~ f) 2 tan $ g ] 


r 


C 


1 - (2f - f 2 ) 

.1 - ( 2 f - f 2 ) cos 2 4>' s 


Compute the height above the reference ellipsoid, 


H s " [r 2 - r 2 sin 2 (0 S - <^)] M - r c cos (<p s - 0 g ) 


A^s 


sin 


-l 



sin (0 S - 0s) , 


7 T 7 T 

~ ~2 — — ~2 


with the declination given by, 

d = 4>§ + A 0g 

and from spherical trigonometry we have for the argument of latitude 
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u 


sin -1 


( sin d \ 
sin i J ' 


thus we have for the true anomaly 


v - u - g 

where the time of the geodetic latitude is determined from Brouwer (1959) p. 395. 
With the time we update the g and the process is continued until; 

/g i+ i - gi/ < e 

where e is some preassigned, small positive number. 


CALLING SEQUENCE 

CALL GTRACE (TAO, TO, REV, GDL, DLONG, HT, TGDL, DGDL, KGDL, 
ADGDL, NSGDL, DRG, NQ, SATID, STAR, TAR) 


INPUT/OUTPUT 

Arguments 


I/O 

Variable 

Description 


TAO 

Tq — start time of reference ellipsoid 


TO 

t 0 — epoch time 


REV 

Number of orbital revolutions 

0 

GDL(360) 

Geodetic latitudes (</) 

o 

DLONG(360) 

Longitude of geodetic meridian east of ascending 
nodal meridian ) 

o 

HT(360) 

Height above meridian of satellites geodetic meridian 

<»y 

o 

TGDL(360) 

Time of geodetic latitude (t^) 

I 

DGDL 

Increment in geodetic latitude by degrees (A 4 >) 

o 

KGDL 

Number of GDL's in the interval ascending node to 
north point 

o 

ADGDL(6) 

ADGDL(l-3) — ty, A.y, hty, 
ADGDL(4-6) - t , \ , ht 
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Arguments (continued) 


I/O 

Variable 

Description 

o 

NSGDL(6) 


NSGDL(l-3) - t N , \ N , ht N , 
NSGDL(4-6) — t s , >v s , ht s 

I 

DRG(l) 

1 


t 0 , t 1# t 19 


DRG(20) 

DRG(21) 

: 

" ^2, O’ ^2,1’ ••••> ^2,19 


DRG(40) 

- 



DRG(41) 

1 



1 

DRG(60) 

. 

‘ ^3,0* ^3,1’ ••••» N 3t 19 

I 

NQ 


Number of drag inputs 

I 

SATID(ll) 


SATID(l) = satellite identification number 

SATID(2) = reference year 

SATID(ll) = day count of reference date 

o 

STAR(300) 


'*» indicates satellites in sunlight, ' * not in sunlight 

o 

TAR(6) 


TAR = BLANK, or TAR = *, indicates if ascending, 




descending nodal crossings, or north and south 
points are in or out of shadow 


Data Statements 


Variable 

Definition 

IT 

IT /10/ max. no. of iterations 

TG 

tg /0.0/ Greenwich epoch time 

ASTRX 

/'*'/ 

BLANK 

/' '/ 
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Common 


I/O 

Block 

Variable 

I 

POOL 

TABLE (9) - 1/f, TABLE(24) = 2, TABLE (31) =t Q l 



TABLE(69) = co e , TABLE(59) - min/cut 

I 

OSCELE 

ORB(6) = a, e, i, g, h, 1 

I 

SECPRM 

DPELE(6) = a", e", i", g", h", 1" 

I 

PERTL 

LDGPT = ldgPt, E 0 ,t 0 

I 

NSMAP 

^minN* ^N> ^EN » » ^N’ ^rainS* * ^ES» ^ S 

I 

DHMNS 

min N , min s , hr N , hr s , Dy , Dy 

I 

DOTELE 

l 0 dot 

I 

NOD 

L0D = 1 q D 

I 

NODMAP 

^Q, t min y» kjj, hty, Ey 

I 

DHMNOD 

min n , min^ hr n , hr y , Dy Q , Dyy 


CALLED BY 
MAIN 

CALLS 

CALL BRWORB 
CALL SSPTHT 
CALL SDFWOE 
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HMSTOR 


Hours, Minutes, Seconds to Radians 


PURPOSE 

To convert hours, minutes and seconds to radians. 


METHOD 

Rad = (sec/60 + minutes)/60 + hours * radians per hour 


CALLING SEQUENCE 

CALL HMSTOR (HH, FM, SS, RAD) 


INPUT/OUTPUT 

Arguments 


I/O 

Variable 

Description 

I 

HH 

Hours 

I 

FM 

Minutes 

I 

SS 

Seconds 

o 

RAD 

Radians 


Common 


I/O 

Block 

Variable 

I 

BPOOL 

TABLE (65) - RADPHR (rad/hr) 


CALLED BY 
MAIN 
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HMSTOR Flowchart 
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INTPLO 


Backward Difference Interpolation 


PURPOSE 

To interpolate for the elements a p , e p , i p , l p , g p , h p when given an observation 
time between two times on the perturbation tape. 


METHOD 


Backward Difference Interpolation 


> times from the perturbation tape 


f(t 6 ) 
f(t 5 ) 
f(t 4 ) 
f (t 3 ) 
f (t 2 ) 
f(t x ) 


elements from the perturbation tape at time t i5 i.e., if we 
> are interpolating for a pert then f(t 6 ) = a 6 , f (t 5 ) = a 5 , 
f (t 4 ) = a 4 , 


t = observation or request time 

At = time increment from perturbation tape 
A K f (t.) = the K th difference of the elements 


A K- 1 f(t. +1 ) - A K " 1 f ( t i ) 
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where 


A° f ( t i ) - f( ti ) 



(elements)^ = f (t 6 ) + Af (t s ) (ttO) + A 2 f (t 4 ) ^ tt0) 2 (ttl > 

+ A 3 f(t 3 ) \ ( ttQ ) (ttl) ( 1 1 2) 1 + a4 f r ( ttO) (ttl) (tt2) (tt3)~ 


, A* V r (ttO) (ttl) (tt2) (tt3) (tt4)l 

U [ 120 

where 


ttO 

= t - t 6 

ttl 

= At + ttO 

tt2 

= At + ttl 

tt3 

= At + tt2 

tt4 

= At + tt3 
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CALLING SEQUENCE 


CALL ENTPLO (TIME, A, B, TSUBO, DELTA) 
INPUT/OUTPUT 


Arguments 


I/O 

Variable 

Description 

I 

TIME 

Observation time 

I 

A(6, 7) 

Array of observation times and elements from the 
perturbation tape 

o 

B(6) 

Interpolated elements from perturbation tape for 
observation time — a p , e p , i p , l p , g p , h p 

I 

TSUBO 

Sixth time in the time element array A 

I 

DELTA 

Time increment between 5 times on the perturbation 
tape 


Definition of Array A (I, J) 


\ 

1 

2 

3 

4 

5 

6 

1 

*1 

t 2 

* 3 

t 4 

^ 5 


2 

a i 

a 2 

a 3 

a 4 

a? 

a 6 

3 

e i 

e 2 

e 3 

e 4 

e 5 

e 6 

4 

i i 

i 2 

* 3 

U 

is 

i 6 

5 

li 

1 2 

1 3 

U 

Is 

1 6 

6 


g 2 

§3 

g 4 

g S 

g6 

7 

hi 

h 2 

h 3 

h 4 

hs 

h 6 


CALLED BY 
PERTFO 
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INTPLO Flowchart 













JDSCUT 


Julian Day — Seconds to Canonical Unit of Time 


PURPOSE 

To convert Julian days (number of days from date of reference to date of obser- 
vation) and seconds to canonical units of time. 


METHOD 

CUT = DAYJ * SECDAY/SECCUT + SS/SECCUT 
CALLING SEQUENCE 


CALL JDSCUT (DAYJ, SS, CUT) 

INPUT/ OUT PUT 
Argument 


I/O 

Variable 

Description 

I 

DAYJ 

Julian days 

I 

SS 

Seconds 

o 

CUT 

Julian days and seconds in canonical units of time 


Common 


I/O 

Block 

Variable 

I 

BPOOL 

TABLE (5) = SECCUT 

I 


TABLE (30) - SECDAY 


CALLED BY 


DRAGLD 

TIMETB 









JDSCUT Flowchart 






JULHMS 


Julian Days — Seconds to Julian Hours, Minutes and Seconds 


PURPOSE 

To convert Julian days (number of days from date of reference to date of obser- 
vation) and seconds to Julian days, hours, minutes and seconds. 


CALLING SEQUENCE 

CALL JULHMS (DAYS, SEC, RF, DAY, HH, FM, SS) 


INPUT/OUTPUT 

Arguments 


I/O 

Variable 

Description 

I 

DAYJ 

Julian days 

I 

SEC 

Seconds 

I 

RF 

Rounding factor (added to SEC) 

o 

DAY 

Julian days 

o 

HH 

Hours 

o 

FM 

Minutes 

0 

SS 

Seconds 


CALLED BY 

SATOR 

TIMETB 
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KEPLRI 

Solution of Kepler's Equation for Eccentric Anomaly 


PURPOSE 

To solve Kepler's equation for eccentric anomaly given mean anomaly and 
eccentricity by the Miles Standish algorithm. 


METHOD 

Kepler's equation for eccentric anomaly M = E + e sin E is solved using the 
Miles Standish algorithm. It is an iterative process dependent upon a tolerance 
value and a maximum number of iterations. Given the mean anomaly, M, and 
the eccentricity, e, the algorithm for computing the eccentric anomaly, E, will 
be: 

1. Set error code = 0 

Set limit of number of iterations, MAX = 10 

2. Set E = 0 

If M = 0, go to Step 13 
If M / 0, go to Step 3 

3. E Q = M + e sin M 

Set number of iterations = 1 


4. 

5. 

6 . 

7. 

8 . 
9. 

10 . 


F = E Q - (e sin E Q ) - M 
D = 1.0 - [e cos (E 0 - 0.5F)] 

E = E 6 - F/D 

If |E 0 - E| - TOL < 0, go to Step 13; otherwise continue to Step 8 
Add 1 to number of iterations 

If (number of iterations - MAX) < 0, continue; otherwise go to Step 12 
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11. Return to Step 4 


12. Set error code = 4 

13. Modulo E by 2 7T 

14. Return to calling program 

The limit of iterations through Steps 4 to 11 is 10. Thus MAX = 10. If this 
number is exceeded, the error code is set to 4. 

TOL is the tolerance at which the last significant digit of the difference between 
the previous calculated eccentric anomaly and the present calculated anomaly is 
allowed. TOL allows an error of ± .05 x 10 " 10 . 


CALLING SEQUENCE 

CALL KEPLRI (MA, ECC, ERRC, E2) 

INPUT/OUTPUT 

Arguments 


I/O 

Variable 

Description 

I 

MA 

Mean anomaly 

I 

ECC 

Eccentricity 

o 

ERRC 

Error code 



= 0, convergence 



/ 0, no convergence 

o 

E2 

Eccentric anomaly 


CALLED BY 

BRWORB 

ELCONO 
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LAGRNO 

Lagrange's 3-Point Interpolation 


PURPOSE 

To interpolate using Lagrange's three-point interpolation. 
METHOD 

Term 1 - YO * (X - XI) * (X - X2) / (XO - XI) * (XO - X2) 

Term 2 = Y1 * (X - XO) * (X - X2) / (XI - XO) * (XI - X2) 

Term 3 = Y2 * (X - XO) * (X - XI) / (X2 - XO) * (X2 - XI) 

Y = Term 1 + Term 2 + Term 3 

CALLING SEQUENCE 


CALL LAGRNO (X, Y, XO, YO, XI, Yl, X2, Y2) 

INPUT/OUTPUT 

Arguments 


I/O 

Variable 

Description 

I 

X 

Input request variable 

o 

Y 

Computed F(X) for output 

I 

XO 

First input variable 

I 

YO 

F(X0) 

I 

XI 

Second input variable 

I 

Yl 

F (XI) 

I 

X2 

Third input variable 

I 

Y2 

F(X2) 


CALLED BY 
SATOR 
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NODALX 


PURPOSE 

To determine nodal crossing times of an earth satellite. 


METHOD 

A two body solution' is determined for the nodal crossing, the Brouwer theory is 
then used to update the now perturbed two-body elements in order to obtain the 
osculating elements corresponding to the two-body solution. 


FORMULATION 

The argument of perigee (g) is the angle between the direction of perigee and 
the ascending node. 

The true anomaly F, is the angle between the direction of perigee and the radius 
vector of the body. 

From these definitions it follows that: 


f. = 

1 


'2tt - g. at the ascending node 
. 77 - g i at the descending node 


and the eccentric anomaly (E) can be obtained from the relations: 


cos E 

1 


cos v + e 

1 + e cos f . ’ 
1 


sin E j - 


Y 1 - e 2 sin f. 
1 + e cos f . 


(1) 


( 2 ) 


where the time of the event is determined from Brouwer [2] p. 395. with this 
time we update the g in Equation (1) and the process is continued until: 

I Bi+l " ill 5 e ( 3 ) 


where e is some preassigned, small positive number. 
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CALLING SEQUENCE 


CALL NODALX (REV, DRG, NQ, TACUT, SATID, TO, IDAD) 

INPUT/OUTPUT 

Arguments 


I/O 

Variable 

Description 

I 

REV 


Number of orbital revolutions from epoch 

I 

DRG(l) 

>l 

1 

V t 0’ t l’ » t 19 


DRG(20) 

DRG(21) 

I 

J 

1 

^2,0’ ^2, 1» • • • •» ^ 2 . 19 


DRG(40) 

DRG(41) 

I 

J 

1 

1 

' ^3,0’ ^3, P ••••* ^ 3,19 


DRG(60) 

J 


I 

NQ 


Number of drag inputs 

o 

TACUT 


Time of ascending nodal crossing from epoch 

I 

SATED(ll) 


SATID(l) = Satellite identification number 

SATID(2) = reference year 

SATID(ll) = day count of reference date 

I 

TO 


Epoch time 

I 

IDAD 


= 0 , time of ascending and descending nodal crossing 




determined 




= 1 , time of ascending node determined 
= 2 , time of descending node determined 


Common 


I/O 

Block 

Variable 

I 

BPOOL 

TABLE (16) = deg/rad 
TABLE (23) = 2 tt 
TABLE (34) = tt 
TABLE (31) = TOL 
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Common (continued) 


I/O 

Block 

Variable 

I 

BPOOL 

TABLE (69) = co e 
TABLE (59) = min/cut 
TABLE (24) = m 2 

I 

OSCELE 

ORB(6) = a, e, i, g, h, 1 

I 

SECPRM 

DPELE(6) = a", e", i", g", h'', 1" 

I 

NOD 

LOD = 1 0 D 

I 

PERTL 

LDGPT = l d „Pt 

o 

NODMAP 

T minf2> ht £2> tmin ht G’ 

o 

DHMNOD 

DATA 

min 0 , min^, hr n , hr^, Dy n , Dyy 
IT = 10 
tg = 0. 0 DO 


CALLED BY 

MAIN 

GTRACE 


CALLS 

CALL BRWORB 
CALL SSPTHT 
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NSPT 


PURPOSE 

Routine to determine north point (maximum satellite geodetic latitude), and 
south point (minimum satellite geodetic latitude). 


METHOD 

A two body solution is determined for the north-south point crossing, the 
Brouwer theory is then used to update the now perturbed two-body elements in 
order to obtain the osculating elements corresponding to the times of the 
north- south point crossing. 


FORMULATION 

Two body analysis leads us to write 

r 7 t 

~2 - gj at the North point 

f ‘ = 3. „ 

~ 2 — g i at the South point 

and the eccentric anomaly (E) can be obtained from the relations: 


cos E 

1 


cos v + e ^ ^ s * n ^ i 

1 + e cos f . ’ S * n 1 + e cos f. 

1 1 


(2) 


where the time of the event is determined from Brouwer (1959) p. 395. With 
this time we update the g in Equation (1) and the process is continued until: 

!g i+ l ■ Bil 1 e 


where e is some preassigned, small positive number. 
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CALLING SEQUENCE 


CALL NSPT (REV, DRG, NQ, SATID, TO, IDNS) 


INPUT/OUTPUT 


Arguments 


I/O 

Variable 

Description 

I 

REV 

Number of orbital revolutions 

I 

DRG(l) 

f ^0 ’ *1* ^19 


DRG(20) 

DRG(21) 

1 

J 

l N 2 _ 0 » N 2> N 2> ig 


DRG(40) 

DRG(41) 

1 

J 

( N 3, O’ N 3, 1 * N 3, 19 


DRG(60) 

J 

I 

NQ 

Number of drag inputs 

I 

SAT ID (11) 

SATID(l) = Satellite identification number 

SATID(2) = reference year 

SATID(ll) = day count of reference date 

I 

TO 

Epoch time 

I 

IDNS 

= 0, time of north point and south point determined 
= 1, time of north point determined 
= 2, time of south point determined 


Common 


I/O 

Block 

Variable 

I 

BPOOL 

TABLE (31) = TOL 



TABLE(69) = w e 



TABLE (59) = min/cut 



TABLE (24) = p 2 

I 

OSCELE 

ORB(6) = a, e, i, g, h, 1 

I 

SECPRM 

DPELE(6) = a”, e", i”, g", h", 1" 

I 

PERTL 

LDGPT = l dg Pt, E 0 , 1 0 

o 

NSMAP 

^minN* ^N» ^EN> ^N» Ht N , t m j nS , X s , X ES , 0 S , Ht s 
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Common (continued) 


I/O 

Block 

Variable 

o 

DHMNS 

DATA 

min N , min s , hr N , hr s , Dy N , Dy s 
IT = 10 
tg = 0.0 DO 


CALLED BY 

MAIN 

GTRACE 


CALLS 

CALL BRWORB 
CALL SSPTHT 
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NSPT Flowchart 














PAGE1 


First Page Print 


PURPOSE 

To write elements, drags, earth constants, and harmonics on tape. 


CALLING SEQUENCE 

CALL PAGE1 (EPOCH, SR NAME, SATID, OUTPUT, NQ, DRAGDT, DRG, 
NTPQ, DATTIM, CDRAG, KDELT, NJ, ELEMO, RUNID) 


INPUT/OUTPUT 

Arguments 


I/O 

Variable 

Description 

I 

EPOCH(IO) 

EPOCH(l-7) — epoch satellite ID and time 
EPOCH(8) = type of the epoch elements 
EPOCH(9) = pass number 
EPOCH (10) = perturbation option 
= 0 no perturbation 
= 1 use perturbation tape 

I 

SRNAME(3) 

Name of satellite 

I 

SATED(IO) 

Reference sat ID and time 

I 

OUTPUT (6) 

Converted epoch elements 

I 

NQ 

Number of drag inputs 

I 

DRAGDT (40) 

Array of packed drag date and time 

I 

DRG(60) ' 

Array of drag time and parameters 

I 

NTPQ 

Number of column times (or cards) 

I 

DATTIM(112) 

Array of packed column date and time 

I 

CDRAG(112) 

Array of column drag parameters 

I 

KDELT(4) 

KDELT(l) and (3) = number of columns to be 

computed per card 

KDELT(2) and (4) = column At in minutes 

I 

NJ 

Number of column times plus one 

I 

ELEM0(6) 

Epoch elements 

I 

RUNID (8) 

Run identification 
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Common 


I/O 

Block 

Variable 

I 

BPOOL 

TABLE (22) = RAD 
TABLE (41) = Mu 
TABLE (42) = FLAT 
TABLE (49) - ROT 
TABLE (58) = J 
TABLE (61) = K2 
TABLE (62) = K3 
TABLE (63) = K4 
TABLE (64) = K5 


CALLED BY 


MAIN 



CO 

cn 




PAG El Flowchart 









CO 

CT> 




PAGE1 Flowchart (continued) 











PERTFO 


Complementary Perturbations 


PURPOSE 

To read the complementary perturbation tape for the Brouwer Orbit Generator. 


CALLING SEQUENCE 

CALL PERTFO (PLN, SATID, TIME, KMULT, B, IERR) 
INPUT/OUTPUT 


Arguments 


I/O 

Variable 

Description 

I 

PLN 

Perturbation tape logical input unit 
> 0 read pert tape on that unit 
<0 do no read pert tape 
= 0 error 

I 

SATED (11) 

SATID(l) = reference satellite ID 
SATED(2) = year of reference 
SATID(3) = day count of reference date 

I 

TIME 

Observation time 

I/O 

KMULT 

K multiplier for AL drag computation 

0 

B(6) 

Array of elements from perturbation tape for 
observation time — a p , e p , i p , l p , g p , h p 

o 

IERR 

Error 

= 0 no error 

= 37 error in reading pert tape 
= 38 wrong satellite ID on pert tape 
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Common 


I/O 

Block 

Variable 

I 

BPOOL 

TABLE (5) - SEC CUT 

I 


TABLE (45) = CUTDAY 


CALLED BY 
BRWORB 


CALLS 

DREFOD 

JDSCUT 

INTPLO 
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PERTFO Flowchart (continued) 








POOL 


Constants Pool 


PURPOSE 

To compute frequently used constants. 


CALLING SEQUENCE 
CALL POOL 


INPUT/OUTPUT 

Common 

I/O Block Variable 

I BPOOL TABLE (1) = meters/ft 

(BLOCK TABLE (2) = km/CUL 

DATA) TABLE (4) = km/A.U. 

TABLE (5) = sec/ CUT 

TABLE (9) = l/flattening coefficient 

TABLE(ll) = u> e (earth rotation in rad/sec) 

TABLE (12) = J 2 

TABLE(13) = J 3 

TABLE(14) = J 4 

TABLE(15) = J 5 

TABLE (16) = deg/ rad 

TABLE (17) = deg 1 obliquity 

TABLE(18) = min >• of 

TABLE (19) = sec J eccliptic 

TABLE(21) = km/mi 

TABLE(22) = radius of earth in CUL 

TABLE (24) = GM = y? 

TABLE (30) = sec/day 
TABLE(32) = sec/hr 
TABLE (33) = deg/hr 
O BPOOL TABLE (41) = y. 

TABLE (42) = Flattening coefficient 
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Common (continued) 


I/O Block Variable 

O BPOOL TABLE(43) = B (polar radius of the earth) 

TABLE (44) = CUL/A.U. 

TABLE (45) = CUT/Day 
TABLE (46) = CUT/hr 

TABLE (47) = Convert CUL/CUT to km/sec 
TABLE (48) = Convert CUL/CUT to km/hr 
TABLE(49) = co e (earth rotation) in rad/CUT 
TABLE (50) = mi/CUL 
TABLE (51) = Convert CUL/CUT to mi/hr 
TABLE (5 2) = e 2 

TABLE (53) = e (eccentricity of earth) 

TABLE (54) = X component of U 2 
TABLE (55) = Y component of U 2 
TABLE (56) = Z component of U 2 

U 2 is an orthogonal unit vector in the ecliptic 
plane expressed in the inertial coordinate system. 
U 2 is perpendicular to U t in the direction of 
positive r. 

TABLE(57) = TAUDOT (mean longitude in rad/CUT) 

TABLE (58) = J 

TABLE(59) = min/CUT 

TABLE(60) = Convert CUL/CUT to m/sec 

TABLE(61) = K2 

TABLE (6 2) = K3 

TABLE (63) = K4 

TABLE (64) = K5 

TABLE (65) = rad/hr 

TABLE (66) = km/ft 

TABLE (67) = Obliquity of ecliptic in rad 

TABLE (68) = KMULT for Drag 

TABLE (69) = a> e (earth rotation) in deg/min 

TABLE (70) = CUT/min 

TABLE(71) = X component of U t 

TABLE (7 2) = Y component of Uj 

TABLE(73) = Z component of U x 

U x is an orthogonal unit vector in the ecliptic 
plane, expressed in the inertial coordinate 


system. U L is directed to the vernal equinox. 
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Common (continued) 


I/O 

Block 

Variable 

o 

BPOOL 

TABLE (74) = Tolerance for R* • U 

TABLE (75) = Tolerance for magnitude of RXU 

TABLE (76)-(80) — Not used 


CALLED BY 
MAIN 

CALLS 

DMSTOR 
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PQUV 


PURPOSE 

Introduction of the orthogonal unit vector set P,Q or U,V into the orbit plane 
coordinate system. 


METHOD 

Aligning both fundamental triads and performing the rotations through Cl, i, oo, 
yields the direction cosines of P,Q; substitution of u for co yields the direction 
cosines for the set U,V. 


FORMULATION 

The direction cosines of P, and Q are given by: 

Pj = cos oo cos Cl - sin oo sin Cl cos i , Q x = - sin oj cos Cl - cos w sin Cl cos i 

P 2 - cos oo sin Cl + sin oo cos Cl cos i , Q 2 = - sin oo sin Cl + cos oo cos Cl cos i 

P 3 - sin 00 sin i , Q 3 = cos 00 sin i , 

and the directional cosines of U, and V are obtained by the substitution of u = 

v + 00 for 00 into the above equations, that is, U = P(i,fi,u) and V = Q(i, D,u) 

CALLING SEQUENCE 

CALL PQUV (I, H, G, NU, P, Q, U, V) 

CALL UV (I, G, H, NU, U, V) 

CALL PQ (I, G, H, P, Q) 
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INPUT/OUTPUT 


Arguments 


I/O 

Variable 

Description 

I 

I 

(i) Orbital Inclination 

I 

G 

(co) Argument of Perigee 

I 

H 

(Cl) Longitude of the Ascending Node 

I 

NU 

(v) True Anomaly 

o 

U(i) 

U(l), U(2), U(3) = U = P(i,0,w) 

0 

V(i) 

V(l), V(2), V(3) - V(i, Cl, oj) 

o 

P(i) 

P, Unit Vector taken as pointing Towards Perifocus 

0 

Q(i) 

Q, Unit Vector in the orbit plane advanced to P by a 
right angle in the direction of increasing true 
anomaly. 


Common 


I/O 

Block 

Variable 

I 

BPOOL 

TABLE (31) = e = .lx 10' 11 


CALLED BY 
BRWORB 

CALLS 

VDOT 


105 
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PR EDS 


Brouwer Prediction Table Bulletin 


PURPOSE 

To calculate precise prediction information for the orbit of a satellite. 


METHOD 

The satellite position can be determined approximately by means of a method 
based on the assumption that the quantities a(t), e(t), I(t), g(t), h(t) and n(t) 
are osculating elements. It is necessary to have a set of mean Brouwer ele- 
ments, the epoch of such elements, a set of times (tj with a constant At, and 
corresponding (N 2 L ) * s whose value may be zero or a non-zero quantity. Other 
associated parameters at epoch time are computed. 

For i = 0, compute 


a(t 0 ) 


1 + 7 2 ( 30 2 - 1 ) 


1 - 


"(to) 



e (t 0 ) 


s i e (t 0 ) 


72 

2e" 


(3 6 2 - 1 ) 


1 - T } 3 


M (t 0 ) = t'o + ilCto) 

4 a (t 0 ) 

K = 3n(t 0 ) 


4 l 1 ~ e (t 0 )) 

3n (t 0 ) 
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For i > 0, compute 


T i,i = t i+ i - t A 


T i.i 


r 2, i 2 


^1 ® 1 (* i + l) ® 1 (*i ) 


C 2 ®l(*i + T 2,i) ' ®l(^i) 


SC, 


- C 2Tlii 2 ] 


11,1 [T l.i T 2,i ( r 2,i - 

M i (t) = M 2i (t - ti ) 2 + M 1>i( t - ti ) + M 0ii 
M o.i = C i + for i = 0 M(tj ) = M 0 

M l.i = K + 8«u.i + 2 .t 1 N 2, j -i( t i - tj.O 

a ( t i + i) = a (t 0 ) - K2N 2 i Tj A + Aa 
e ( t i + i) = e(t 0 ) + Sj e (t i + 1 ) - e (t 0 ) - W2N 2i Tj . + Ae 


II 

+ 

HH 

I" + 5 ! I (t 

gCti+i) 

= r'cw 

h (t i+ l) 

= h'(t i + 1 ; 

1 

“ 3 Cc> (t Q ) 

^( t i + l - ti ) 
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Q(t 


1 


i) = 2(t +i _ t ~ -3fi(t 0 ) + 40(1,) - 0(t 2 ) 


Period, nodal (t 0 ) 


2 77 


M l,^> + W (*<*>) 


2 77 

Period, anomalistic (t, ) - "Tj 

•l.i 


-M2,i P a,i 2 


a , i 


Height of perigee - a (1 - e) 1 

Height of apogee = a (1 + e) - 1 


Velocity of perigee 



Velocity of apogee 



t 2 + [2 (2 77 - M 0 2 )] 

TP =. 1 ; 

M, 2 jl + ]/l + [4M 2i2 (277 - M 0i2 )] /M 1>2 2 ) 


From BRWORB, compute 

»i at (t, + T 2?i ) 

S , e, s x I, C,, g\ h\ Aa, A e at (t. + 1 ) 
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CALLING SEQUENCE 


CALL PREDS (TO, TPQ, NJ, CDRAG , DRG, EL, PA, PADOT, PNA, TP, 
NN, MI, ELEMO, NQ, SATID) 


INPUT/OUTPUT 

Arguments 


I/O 

Variable 

Description 

I 

TO 

Epoch time 

I 

TPQ(56) 

Column elements times 

I 

NJ 

Number of columns to be computed 

I 

CDRAG(112) 

Column drag table 

I 

DRG(60) 

Brouwer drag time and parameters table 

o 

EL(8, 56) 

Array of Bulletin prediction elements — a . , e . , I . , 

» M 0 , i ’ M i , i » M 2 , i where i = 1-56 
Array of Bulletin prediction anomalistic periods 

o 

PA (56) 

o 

PADOT (5 6) 

Array of prediction anomalistic period derivatives 

o 

PNA(8) 

PNA(l) = nodal period 
PNA(2) = prediction « 

PNA(3) = prediction D 

PNA(4) = perigee height 

PNA(5) = apogee height 

PNA(6) = geocentric latitude of perigee 

PNA(7) = perigee velocity 

PNA (8) = apogee velocity 

o 

TP 

Time of perigee 

o 

NN 

Index for column times 

0 

MI 

Index for column prediction elements 

I 

NQ 

Number of drag inputs 

I 

SATED(ll) 

SATED (1) = reference satellite ID 

SATED(2) = year of reference 

SATED (11) = day count of reference date 
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Common 


I/O 

Block 

Description 

I 

BPOOL 

TABLE (24) = BK 
TABLE (34) = tt 
TABLE (41) - p. 

I 

ETAP 

ETA3 - V 3 ) 
ETA4 -f? 4 ) 
ETA6 - V 6 ) 

I 

GMPR 

GM2 - 7 2 ) 

I 

LPPRM 

DELIE- S,e 
DEL1I - Sji 
L. - Lj 
GP - g 1 
HP - h‘ 

I 

SECPRM 

DPELE(6) - Brouwer mean elements 

I 

DOTELE 

LODOT - L 0 

I 

THETAP 

M1P3T2 - 3<9 2 - 1 

I 

DELKEP 

DKEP(l) - Aa pert 
DKEP(2) - Ae pert 

I 

PRTKEP 

PKEP(3) - C 0 +A C 


TAPE OUTPUT 
Prediction Space Elements 


Epoch: Calendar Date 
Epoch: Julian Date for Space 



^ 2 

t 3 

Period 

minutes 

P 1 

P 2 

P 3 

Period derivative 

microday s/day 

?! 

P 2 

P 3 

Eccentricity 


e i 

e 2 

e 3 

Inclination 

degrees 

*1 


h 

Right ascension of ascending node 

degrees 

Q x 

Q 2 

Q 3 

Argument of perigee 

degrees 

"i 

o) 2 

M . 3 

Mean anomaly 

degrees 

M ! 

m 2 

M 2 

Semi-major axis 

earth radii 

a i 

a 2 

a 3 
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CALLED BY 


MAIN 

CALLS 

BRWORB 

REDUCE 





























PRINT 


Brouwer Prediction Print 


PURPOSE 

To write pertinent prediction information for the orbit of a satellite. 


CALLING SEQUENCE 

CALL PRINT (NC, NQ, SRNAME, EPOCH, ETIME, NN, MI, PA, PADOT, 
EL, PNA, TPQ, SATID, SECDRG, DATTIM, DRG, OSCO, 
PV, ISSUE, NEPASS, TO, ELEMO) 


INPUT/OUTPUT 

Arguments 


I/O 

Variable 

Description 

I 

NC 

Number of columns to be printed 

I 

NQ 

Number of drag inputs 

I 

SRNAME (3) 

Satellite name 

I 

EPOCH(IO) 

Epoch time and information 

I 

ETIME 

Epoch hours, minutes, seconds in seconds 

I 

NN 

Index for column times 

I 

MI 

Index for column prediction elements 

I 

PA(56) 

Array of anomalistic periods 

I 

PADOT(56) 

Array of anomalistic period derivatives 

I 

EL(8, 56) 

Array of Bulletin prediction elements — a. , e. , I. , 
Q., co., M n . , M , , , M, . where i = 1, 56 
PNA(l) = nodal period 
PNA(2) = prediction h 
PNA(3) = prediction co 
PNA(4) = perigee height 
PNA (5) = apogee height 
PNA(6) = geocentric latitude of perigee 
PNA(7) = perigee velocity 
PNA(8) = apogee velocity 

I 

PNA(8) 

I 

TPQ(56) 

Array of column elements times 
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Arguments (continued) 


Vo 

Variable 

Description 

I 

SATED(ll) 

SATID(l) = reference satellite ID 
SATID(2) = year of reference 
SATID(ll) = day count of reference date 

I 

SECDRG(112) 

SECDRG(l), (3), . . (Ill) - seconds from columns 
time 

SECDRG(2), (4), . . (112) = total column time 
(hours, minutes, seconds) in seconds 

I 

DATTIM(112) 

Packed column date and time 

I 

DRG(60) 

Drag time and parameters table 

I 

OSCO(6) 

Epoch Brouwer osculating elements 

I 

PV(6) 

Epoch Brouwer position and velocity vectors 

I 

ISSUE (4) 

Date of issue 

I 

NBPASS 

Initial pass number 

I 

TO 

Epoch in Canonical Unit of Time 

I 

ELEM0(6) 

Epoch elements — a, e, i, g, h, m 


Common 


I/O 

Block 

Variable 

I 

BPOOL 

TABLE(l) =meters/ft 
TABLE (2) = km/CUL 
TABLE (5) = sec/ CUT 
TABLE (16) = deg/rad 
TABLE (23) = 2 tt 
TABLE( 25) = knots/mi/hr 
TABLE (26) = km/N.M. 

TABLE (28) = GM (km 3 /sec 2 ) 

TABLE (30) = sec/day 
TABLE (34) =tt 
TABLE (35) = min/day 
TABLE (41) = fi 
TABLE (45) = cut/ day 

TABLE(49) = u e (earth rotation) in rad/ CUT 
TABLE (50) = mi/CUL 
TABLE (51) = CUL/CUT to mi/hr 
TABLE (57) = TAUDOT 
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Common (continued) 


I/O 

Block 

Variable 

I 

BPOOL 

TABLE (59) = min/CUT 
TABLE (60) = CUL/CUT to m/sec 
TABLE (65) = rad/hr 
TABLE (66) - km/ft 

TABLE (67) = obliquity of ecliptic in rad 

I 

RADIAN 

TAU = t, satellite degrees, minutes, seconds in rad 
AMBDA = k, satellite. hour, minutes, seconds in rad 


CALLED BY 
MAIN 


CALLS 

DREFOD 

REDUCE 

DATANO 

SPACEL 
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PRINT Flowchart 


















REDUCE 


Reduce Angle 


PURPOSE 

To reduce an angle between 0 and 2 tt. 


CALLING SEQUENCE 
REDUCE (Z, PI2) 

Note that REDUCE is a function. 


INPUT/OUTPUT 

Arguments 


I/O 

Variable 

Description 

I 

Z 

Angle to be reduced 

I 

PI2 

2 rr - 360° in radians 

o 

REDUCE 

Reduced-angle Z 


CALLED BY 

BRWORB 

DAF 

DRAG 

PREDS 

PRINT 

SATOR 
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C 


ENTER \ . 



REDUCE = 

REDUCE ) 1 

Y = 1 2 | 


Y Modulo 2 7T 



Negative 


^ RETURN 


REDUCE Flowchart 








SDFWOE 


Sunlight Determination 


PURPOSE 

To determine whether a satellite is in sunlight or darkness (due to the earth's 
shadow) at a given time. An option is available to consider the effects of an 
oblate earth in making this determination. 


METHOD 

Given: 


f 

TIME 


= the position vector of the satellite in the inertial coordinate 
system, measured in CUL, and having components X, Y, 
and Z. 

= the longitude of the sun on reference date, in radians. 

= the motion of tau, in radians per CUT. 

= t - t 0 = the time in CUT, measured from reference date, at 
which the sunlight determination is to be made. 


Uj and U 2 = orthogonal unit vectors in the ecliptic plane, expressed in the 
inertial coordinate system. U j is directed to the vernal 
equinox and U 2 is perpendicular to U t in the direction of 
positive r (U x = 1, 0, 0; U 2 = 0, cos e, sin e where e = 
obliquity of ecliptic). 

f = the flattening coefficient of the ellipsoid of reference. 

Compute: 


r* = the unit satellite position vector 
T = r + f (t - tjj) 

U = Uj cos T + U 2 sin T, having coordinates u, v, w 
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If | r x u | < Tj, where T t = 1 - f (Z + W y'r 2 - l) 2 , or a constant, and if 
r* • U < T 2 , then the satellite is in darkness. Otherwise, it is in sunlight. 

CALLING SEQUENCE 

CALL SDFWOE (TIME, R, IX, RMAG) 


INPUT/OUTPUT 

Arguments 


I/O 

Variable 

Description 

I 

TIME 

Time at which sunlight determination is made 

I 

R(3) 

Position vector of the satellite — X, Y, Z 

o 

IX 

Sunlight Determination 



= 0, satellite is in darkness 



= 1, satellite is in sunlight 

I 

RMAG 

Magnitude of position vector R 


Common 


I/O 

Block 

Variable 

I 

BPOOL 

TABLE (42) = F 

I 


TABLE(54) = X component of U 2 

I 


TABLE(55) = Y component of U 2 

I 


TABLE (56) = Z component of U 2 

I 


TABLE (57) = TAUDOT (f) 

I 


TABLE (71) = X component of U x 

I 


TABLE (72) = Y component of Uj 

I 


TABLE(73) = Z component of Uj 

I 


TABLE (74) = T2 (tolerance for R* • U) 

I/O 


TABLE(75) = T1 (tolerance for magnitude of RXU) 

I 

RADIAN 

TAU (r) 
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CALLED BY 


GTRACE 


CALLS 

VECTOR 
VADD 
VDOT 
VMAG 
VPROD 
VUNIT 
VC ROSS 
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SDFWOE Flowchart 


















SPACEL 


Spacel Bulletin 


PURPOSE 

To provide key space element information in a condensed form. 


CALLING SEQUENCE 

CALL SPACEL (ISSUE, SPAC, DRG, NQ, TO, TPQ, NBPASS, NSPACE, 
SATID) 


INPUT/OUTPUT 


Arguments 


I/O 

Variable 

Description 

I 

ISSUE (4) 

Issue Date 

I 

SPAC (487) 

May consist of 54 of the following set of elements 
depending on the number of elements columns to 
be computed 
SPAC(l) = time 
SPAC (2) = anomalistic period 
SPAC (3) = derivative of anomalistic period 
SPAC (4) = perigee height 
SPAC (5) = apogee height 
SPAC (6) = inclination 

SPAC (7) = right ascension of ascending node 
SPAC (8) = argument of perigee 
SPAC (9) = mean anomaly 

Last SPAC(N) = 9999999999999999 (end of spacel) 

I 

DRG(60) 

Drag time and parameters table 

I 

NQ 

Number of drag inputs 

I 

TO 

Epoch time 

I 

TPQ(56) 

Table of column elements times 

I 

NBPASS 

Request pass number 

I 

NSPACE 

Number of days from space reference date (Sept. 18, 
1957) to issue date 
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Arguments (continued) 


; I/O 

Variable 

Description 

I 

SATID(ll) 

SATID(l) - satellite ID number 
SATED(2) - year of reference 
SATID(ll) - day count of reference date 


Common 


I/O 

Block 

Variable 

I 

OSCELE 

OS(6) = Brouwer osculating elements 

I 

BPOOL 

TABLE (2) = Km/CUL 
TABLE(16) = deg/rad 
TABLE (23) - 2 tt 
TABLE (24) = BK 
TABLE (34) = tt 
TABLE (40) = NSPAC 

= 0, no spacel osculating epoch elements output 

= 1, spacel osculating epoch elements output 
TABLE (41) = fj. 

TABLE(59) = min/CUT 


TAPE OUTPUT 


GOOOARD SPACE FLIGHT CENTER 
SPACEL BLLLETIN 
MAR 3 » 1 9 71 

j.d.s. asis 


SATCESC NAMF SL NRDN SBDAT 
.1.0. S. PERAN MOO PER CERI 
»JT2W 10,000 MIN MICROD/D 

mno xuccc.ccc xic.coc 


INJ.O. 5. R VNEI OMAG 
PHTEOR APHTEOR INMCD 
KMXlO XMX10 1000G 
2Se.25 <378166 X1COO 


KG/ M2 RFDxwCM FCW CCNC 
RAANOD , A KG P E P MANCM C 
DEG PS CtGRS CEGRS S 
X1000 X lo 0 3 XIOOO C 


68C6602 INJUN5 UW 333E ASIA 3S1A 
A <50600 Ollf 2E8eS3-COCOCAC5 CC67C1 
A913C0 C11E2E8A8A— CCCOC<28 CC67C3 
A92C<T Cll E288CA2-CCCCCA67 CC67C9 
A9C600 011816121 1-CCC0C393 0C6618 


397 112921292 

C025327 8C668 3A5393 092867 207521A 
CC25325 80668 3A0103 07b81A 2850103 
CC2531 8 8C668 33A81A 064753 0026103 
00253A 8C 8C665 S3A5 3 9 1 COS j 1 27 J 20 7237 A 
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The first line of data contains the satellite identification number, satellite name 
and other information from the Spacel input card. 

The following lines of data each contain: 

1. Date measured from the Julian Date of Space 

2. Anomalistic period in minutes, modulo 10,000, X10 6 

3. Period derivative, microdays/ day, X10 4 

4. Perigee height relative to equatorial radius, kilometers, X10. 
Flattening = 1/298.25 

5. Apogee height relative to equatorial radius, kilometers, X10. 
Equatorial radius = 6378.166 

6. Inclination in degrees, modulo 100, X1000 

7. Right ascension of ascending node in degrees, X1000 

8. Argument of perigee in degrees, X1000 

9. Mean anomaly in degrees, X1000 

10. Check sum of line (last digit), modulo 10. 

The last test line is computed by using the osculating elements at the start 
request time. 


CALLED BY 
PRINT 


CALLS 

CKSUM 

ZERO 

BRWORB 
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SPACEL Flowchart 
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SPACEL Flowchart (continued) 



SSPTHT 


PURPOSE 

This routine determines sub-satellite point and height, i.e., a. transformation 
from cartesian coordinates to position in the latitude- longitude coordinate; 
system . 


METHOD 

As a first approximation the geocentric latitude is set equal to the satellite's 
declination, an iterative procedure is then employed to determine the geodetic 
latitude and height. 


FORMULATION 


Sin S 


Z 
r ’ 


7 7 



< 


77 

~2 


k' E - a - e g - W e t , \ E 


mod (\ E , 2 77), 0 £ < 2n 


Set = S, where S has already been determined, and continue calculating with 


a. 


1 - ( 2 f - f 2 ) 

1 - (2f - f 2 ) cos 2 0, 


( 1 ) 


0 S - tan" 


tan 0c 


(1 - f) 5 


77 77 

" ~2 - - ~2 


H s = [r 2 - r c 2 sin 2 (0 - 0s)] 1,4 “ r c cos (0 S - 0g) 


H c 


A0; 


sin 


-l 


sin (0 S - 0g) 


77 77 

~2 — ^ — ~2 
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Recalculate 4> g = § - A cp' s and return to Eq. (1) . Repeat this loop until no 
longer varies. This process is exact and rapidly convergent, and at the same 
time yields the ground trace geodetic latitude, 4>. 


CALLING SEQUENCE 

CALL SSPTHT (WEDT, R, RMAG, RTASC, DEC, GEODL, GEOCL, LONG, 
ELONG, HEIGHT) 


INPUT/OUTPUT 

Arguments 


Vo 

Variable 

Description 

i 

WEDT 

oj e (t - t 0 ) Earth's Rotation from Start Date 

i 

R(3) 

X, Y, Z Satellite Position Vector 

i 

RMAG 

Magnitude of Sat. Radial Vector 

o 

RTASC 

Right Ascension — Deg. 

o 

DEC 

Declination — Deg. 

o 

GEODL 

Geodetic Latitude — Deg. 

o 

GEOCL 

Geocentric Latitude — Deg. 

o 

LONG 

Longitude - -n < LONG < tt measured from the 
Greenwich Meridian 

o 

ELONG 

East Longitude (k E ) — Deg. 0, 2 v measured 
Eastward from Greenwich Meridian 

o 

HEIGHT 

HEIGHT — Above the reference ellipsoid in km 


Common 


Vo 

Block 

Variable 

I 

I 

SDT 

BPOOL 

GST — Greenwich Sidereal Time 
TABLE (31) - e = .1 x 10" 11 
TABLE (9) = 1/f 
TABLE (2) = km/E.R. 

TABLE (42) = f 
TABLE (52) = 2f - f 2 
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CALLED BY 


G TRACE 
NODALX 
NSPT 


CALLS 


None 




SSPTHT Flowchart 











TIMETB 


Brouwer Column Elements Time Table and Drag Load 


PURPOSE 

To handle the input of the column elements times and drags. 


CALLING SEQUENCE 

CALL TIMETB (SATED, NTPQ, KDELT, DATTIM, SECDRG, TPQ, CDRAG, 
NERROR, NJ) 


INPUT/OUTPUT 


Arguments 


I/O 

Variable 

Description 

I 

SAT ID (11) 

SATID(l) = satellite ID number 
SATID(2) = year of reference 
SATID(ll) = day count of reference date 

o 

NTPQ 

Number of column times (or cards) 

o 

KDELT(4) 

KDELT(l) and (3) - Number of columns to be 
computed per card (maximum of 2 cards) 
KDELT (2) and (4) - Column AT in minutes x 100 

o 

DATTIM(112) 

DATTIM(1, 3, 5 . . . Ill) - Packed column date 
DATTIM(2, 4, 6 . . . 112) - Packed column time 

o 

SECDRG(112) 

SECDRG(1, 3, 5 . . . Ill) - Column seconds 
SECDRG(2, 4, 6 . . . 112) - Column hr, min, sec in 
seconds 

o 

TPQ (5 6) 

Column times in CUT 

o 

CDRAG(112) 

Column drag parameters (N 2 q , N 3 q ) 

o 

NERROR 

Error indicator 
= 0 no error 

= - 1 error on column card 

o 

NJ 

Number of column times plus one 
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Common 


I/O 

Block 

Variable 

I 

BPOOL 

TABLE (5) = sec/CUT 
TABLE (30) = sec/day 


CALLED BY 
MAIN 


CALLS 

DREFOD 

JDSCUT 

JULHMS 

DATE 
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UTGST 


PURPOSE 

This routine calculates the Greenwich sidereal time 


METHOD 

Calculate the Greenwich sidereal time at 0 hr U.T. of Date and add the 
amount of rotation in the hours, minutes, and seconds elapsed since 0 hr U.T. 
of Date. 


FORMULATION 

The Greenwich sidereal time at 0 hr U.T. is given by 

0 g 0 = 99?6909833 + 36000?7689 T u + 0?00038708 T u 2 

where the time is measured in centurie as 

J.D. - 2415020.0 
u 36525 

and we have for the Greenwich sidereal time 

(9g = 0g o + w e (t - t^) 

where w e is the earth's rotation in deg/min 

(t - t Q ) is the number of minutes elapsed since 0 hr U.T. 

CALLING SEQUENCE 

CALL UTGST (YY, MM, DY, HR, MIN, SEC) 
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INPUT/OUTPUT 


Arguments 


I/O 

Variable 

Description 

I 


Year 

I 


Month 

I 


Day 

I 


Hour 

I 

MIN 

Minute 

I 

SEC 

Second 


Common 


I/O 

Block 

Variable 

0 

SDT 

GST, GSTO - (0g,'0g o ) 

I 

BPOOL 

TABLE (69) - w e 


CALLED BY 
MAIN 


CALLS 

None 
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UTGST Flowchart 


UVIJK 


PURPOSE 

Maps the radius vector to the satellite from the orbital (orthogonal set u, v , w) 
to the equatorial coordinate system. 


METHOD 

Call subroutine UV (I, G, H, NU, U, V) to obtain u, v unit vectors in the orbital 
plane then calculate the positional vector and the velocity vector in the inertia 
system. 


FORMULATION 

r = ru, r = a(l-e cos E) 
v r = VmI /r e Sin E, v T = 3 /EL ( ,L~ e D 
r = v r u + v T v 


CALLING SEQUENCE 

CALL UVIJK (A, ECC, I, G, H, F, E, R, DR, U, V, RMAG, DRMAG) 

INPUT/OUTPUT 

Arguments 


I/O 

Variable 

Description 

I 

A 

Semimajor axis (e.r.) 

I 

ECC 

Eccentricity 

I 

I 

Orbital inclination 

I 

G 

Argument of perigee 

I 

H 

Longitude of the ascending node 

I 

F 

True anomaly 

I 

E 

Eccentric anomaly 
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Arguments (continued) 


I/O 

Variable 

Description 

o 

R 

Sat. radial vector 

o 

DR 

Sat. vel. vector 

o 

U 

P(I, G & F, H) = r/r 

o 

V 

Q(I, G & F, H) 

o 

RMAG 

r magnitude of radial vector 

o 

DRMAG 

v magnitude of vel. vector 


Common 


I/O 

Block 

Variable 

I 

BPOOL 

TABLE (31) = e = .1 x 10’ 11 

I 

BPOOL 

TABLE (41) = ju 

O 

VTVR 

VT(3) Transverse velocity component vector 

O 

VTVR 

VR(3) Radial velocity component vector 


CALLED BY 
BRWORB 


CALLS 

VDOT 
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UVIJK Flowchart 









VECTOR 


Vector Operations Package 


PURPOSE 

VECTOR performs the following vector operations: 

• Performs addition of two vectors . 

• Computes the cross product of two vectors. 

• Computes the dot product of two vectors. 

• Computes the magnitude of a vector. 

• Computes the product of a scalar and a vector. 

• Performs subtraction of two vectors. 

• Computes the unit vector along a given vector. 

METHOD 

VADD — Performs addition of two vectors. 

The vector sum, S = (S v S 2 , S 3 ), of the two given vectors, A = (A v A 2 , A 3 ) 
and B = (B 15 B 2 , B 3 ), is computed as follows: 

Sx = \ + 

^2 - ^2 + 

s 3 - a 3 + b 3 

VCROSS — Computes the cross product of two vectors. 

The cross product vector, C = (Cj, C 2 , C 3 ), of the two given vectors, A = 
(Aj, A 2 , A 3 ) and B = (Bj, B 2 , B 3 ), is computed as follows: 
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a 2 b 3 a 3 b 2 


Cl = 

c 2 _ A 3 Bj - AjB 3 
c 3 ~ AjB 2 A 2 Bj 


VDOT — Computes the dot product of two vectors. 

The dot product, A • B, of the two given vectors, A = (A x , A 2 , A 3 ) and B = 
(B p B 2 , B3), is given "by: 

A " B - AjBj + A 2 B 2 + A 3 B 3 
VMAG — Computes the magnitude of a vector. 

The magnitude A of the given vector, A = (A v A 2 , A 3 ), is computed as follows: 

A = ( Al 2 + A 2 2 + A /) 54 

VPROD — Computes the product of a scalar and a vector. 

The product vector, P = (P 3 , P 2 , P 3 ), of the scalar, a, and the vector, A = 

(A 3 , A 2 , A 3 ), is computed as follows: 

p i = aA i 

P 2 - a A 2 

P 3 - a A 3 


VSUB — Performs subtraction of two vectors. 

The vector difference, D = (D t , D 2 , D 3 ), of the two given vectors, A = 
(Aj, A 2 , A 3 ) and B = (B p B 2 , B 3 ), is computed as follows: 
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Di - A. l Bj 

D 2 ~ A 2 _ b 2 

d 3 - a 3 - b 3 


VUNIT — Computes the unit vector along a given vector. 


The unit vector A* = (A*, A*, A*) along the vector A = (A 1# A 2 , A 3 ) is given 
by: 


where 


A** 


h. 

|A| 


A* 

a 2 


A 2 

III 


A* 


A 


3 


I A 


I A | = (A 2 + A| + A 3 2 )* 


CALLING SEQUENCE 

CALL VADD (A, B, APLUSB) 

CALL VCROSS (A, B, ACRSSB) 
CALL VDOT (A, B, ADOTB) 

CALL VMAG (A, AMAG) 

CALL VPROD (A, SCALAR, PROD) 
CALL VSUB (A, B, AMNSB) 

CALL VUNIT (A, AUNIT) 
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INPUT/OUTPUT 


Arguments 


I/O 

Variable 

Description 

I 

A(3) 

Input vector on which a vector operation is to be 
performed 

I 

B(3) 

Input vector on which vector operation is to be 
performed 

0 

APLUSB 

Vector sum of A and B 

0 

ACRSSB 

Cross product of vectors A and B 

o 

ADOTB 

Dot product of vectors A and B 

0 

AMAG 

Magnitude of A 

I 

SCALAR 

Scalar input to be multiplied by the vector A 

0 

PROD 

Product of the scalar and vector A 

o 

AMNSB 

Vector difference of A and B 

0 

AUNIT 

Unit vector along vector A 


CALLED BY 

ELCONO 

PQIJK 

PQUV 

SDFWOE 

UVIJK 
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VECTOR Flowchart 












WMAPLD 


World Map Request Load 


PURPOSE 

To handle input for the world map (equator crossings and one orbit ephemeris). 


CALLING SEQUENCE 

CALL WMAPLD (REQUST, START, END, SATID, JDREQ, RUNID, NERROR) 

INPUT/OUTPUT 

Arguments 


I/O 


Variable 


Description 


O 


O 

0 

1 

O 

O 


REQUST(15) 


START 

END 

SATED(ll) 

JDREQ 

NERROR 


REQUST(l) = Year 
REQUST(2) = Month 
REQUST (3) = Day 
REQUST (4) = Hour 
REQUST(5) = Minutes 
REQUST (6) = Seconds 
REQUST(7) = Year 
REQUST (8) = Month 
REQUST (9) = Day 
REQUST (10) = Hour 
REQUST (11) = Minutes 
REQUST(12) = Seconds ^ 

REQUST(13) = Latitude increment 

REQUST(14) = Inclination 

REQUST (15) = Pass number 

Start time of request in CUT 

End time of request in CUT 

SATID(2) = Reference year 

SATID(ll) = Day count of reference date 

Number of days from date of reference to start 

date of request 
Error indicator 

= 0 no error 

= - 1 latitude increment > inclination 


Start Time of 
WMAP Request 


End Time of 
WMAP Request 
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CALLED BY 
MAIN 

CALLS 

DREFOD 

JDSCUT 



149 



WMAPLD Flowchart 






ZERO 


Zero Print 


PURPOSE 

To print leading zeroes of an integer on the IBM 360 computer. 


CALLING SEQUENCE 

CALL ZERO (N, IIN, AREA) 

INPUT/OUTPUT 

Arguments 


I/O 

Variable 

Description 

I 

N 

Number of output digits (including the sign) that is 
desired — maximum is 8 

I 

IIN 

Input number 

o 

AREA(8) 

AREA(l) = Sign of the output number 
AREA(2)-(8) = digits of output number including 
leading zeroes. 


CALLED BY 
SPACEL 
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ZERO Flowchart 






















BLOCK DATA 


PURPOSE 

To assign values to frequently used constants which do not have to be computed. 


OUTPUT 

Common 


Block 

Variable 

BPOOL 

TABLE (1) y meter/ ft 
TABLE (2) - km/CUL 
TABLE (3) = mile/nautical mile 
TABLE (4) = km/A.U. 

TABLE (5) = sec/ CUT 

TABLE(6) = velocity of light (km/sec) 

TABLE (7) = sun mass/ earth mass 

TABLE(8) = earth mass/ moon mass 

TABLE (9) = 1/f (flattening coefficient) 

TABLE(IO) = Pressure of sunlight 

TABLE (11) = co e - earth rotation (rad/sec) 

TABLE (12) = J 2 

TABLE(13) = J 3 

TABLE (14) = J 4 

TABLE (15) = J 5 

TABLE (16) = deg/rad 

TABLE(17) = Deg 

TABLE (18) = Min Obliquity of Ecliptic 

TABLE(19) = Sec 

TABLE(20) = Mean long, of sun (deg/ day) 

TABLE(21) = km/mile 

TABLE (22) = Radius of earth in CUL 

TABLE(23) - 2 tt 

TABLE (24) = GM = p. 2 

TABLE (25) = knots/mi/hr 

TABLE(26) = km/n.m. 

TABLE (27) = Lunar distance (CUL) 
TABLE(28) = GM (km 3 /sec 2 ) 
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Common (continued) 


Block 

Variable 

BPOOL 

TABLE (29) = Solar distance (CUL) 
TABLE(30) = sec/day 

TABLE (31) = Tolerance for orbit generator 
TABLE (32) - sec/hr 
TABLE (33) = deg/hr 
TABLE (34) = -n 
TABLE (35) - min/day 
TABLE (36) = ID for sator code 
TABLE (37) = Radio frequency 
TABLE (38) = Radio frequency 
TABLE (39) = Radio frequency 
TABLE (40) = SPACEL osculating epoch 
elements output option 
= 0 do not compute SPACEL osculating 
elements 

= 1 compute SPACEL osculating elements 





C. COMMON BLOCK VARIABLE DESCRIPTION 


The following section contains the COMMON block descriptions of the COMMON 
areas used in the Goddard Brouwer Orbit Bulletin program. These descriptions 
include the variables contained in these areas, their meaning, and the subroutine 
which defines each variable. 


1. BULLETIN COMMON Blocks 

This section contains a description of all the COMMON areas used in the 
BULLETIN program. 


/BPOOL/ 

COMMON/ BPOOL/ TABLE (80) 


Variable 

Value 

Description 

Program 
Where Defined 

TABLE 


Table of Constants 


TABLE(l) 

.30480 

meters/foot 

BLOCK DATA 

TABLE(2) 

.6378.1660 

kilometers/ Canonical Unit 

BLOCK DATA 



of Length 



TABLE (3) 

1852.0 

meters/nautical mile 

BLOCK DATA 

TABLE (4) 

149598600.0 

kilometers/ astronomical unit 

BLOCK DATA 

TABLE (5) 

806.812418099482 

seconds/ Canonical Unit of Time 

BLOCK DATA 

TABLE (6) 

299792.5 

velocity of light in km/ sec 

BLOCK DATA 

TABLE (7) 

332948.55 

ratio of sun mass/ earth mass 

BLOCK DATA 

TABLE (8) 

81.3 

ratio of earth mass/moon mass 

BLOCK DATA 

TABLE (9) 

298.250 

l/flattening coefficient (f) 

BLOCK DATA 

TABLE(IO) 

.000045 

pressure of sunlight 

BLOCK DATA 

TABLE (11) 

.0000729211510 

w e - earth rotation in rad/ sec 

BLOCK DATA 

TABLE(12) 

.001082480 

J 2 


BLOCK DATA 

TABLE(13) 

- .00000256 

J 3 


BLOCK DATA 

TABLE (14) 

- .00000184 



BLOCK DATA 

TABLE(15) 

- .00000006 

J 5 


BLOCK DATA 

TABLE(16) 

57.29577951308232 

degrees/ radian 

BLOCK DATA 

TABLE (17) 

23.0 

degrees 

Obliquity 

BLOCK DATA 

TABLE(18) 

26.0 

minutes 

> of 

BLOCK DATA 

TABLE(19) 

34.795 

seconds ^ 

Ecliptic 

BLOCK DATA 

TABLE(20) 

.9856470 

mean longitude of sun in deg/day 

BLOCK DATA 

TABLE(21) 

1.609344 

kilometers/ mile 

BLOCK DATA 

TABLE (22) 

1.0 

R - radius of the earth in CULje.r.) 

BLOCK DATA 

TABLE(23) 

6.283185307179586 

2 77 (360° in radians) 

BLOCK DATA 
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/BPOOL/ (continued) 


Variable 

Value 

Description 

Program 
Where Defined 

TABLE(24) 

1.0 

GM = p 2 (Gaussian constant) 

BLOCK DATA 

TABLE(25) 

1.152 

knots/ mile/hr 

BLOCK DATA 

TABLE(26) 

1.852 

kilometers/ nautical mile 

BLOCK DATA 

TABLE (27) 

60.266011 

lunar distance in CUL 

BLOCK DATA 

TABLE(28) 

398603.2 

GM* (km 3 / sec 2 ) 

BLOCK DATA 

TABLE(29) 

234658.04 

solar distance in CUL 

BLOCK DATA 

TABLE(30) 

86400.0 

seconds/ day 

BLOCK DATA 

TABLE(31) 

.000000000001 

tolerance 

BLOCK DATA 

TABLE (32) 

3600.0 

seconds/ hour 

BLOCK DATA 

TABLE (33) 

15.0 

degree s/hour 

BLOCK DATA 

TABLE (34) 

3.141592653589793 

77 (180° in radians) 

BLOCK DATA 

TABLE(35) 

1440.0 

minutes/ day 

BLOCK DATA 

TABLE (36) 

0.0 

indicator for sater code 
= 0 compute sater code 
> 0 no sator code 

BLOCK DATA 

TABLE(37) 

0.0 

radio frequency 

BLOCK DATA 

TABLE (38) 

0.0 

radio frequency 

BLOCK DATA 

TABLE (39) 

0.0 

radio frequency 

BLOCK DATA 

TABLE (40) 

1.0 

indicator for SPACEL osculating 
elements 

= 1 compute osculating elements 
= 0 no SPACEL osculating 
elements 

BLOCK DATA 

TABLE(41) 

1.0 


POOL 

TABLE (42) 

f 

flattening coefficient 

POOL 

TABLE (43) 

(1 - f) * R 

B - polar radius of the earth 

POOL 

TABLE(44) 

(km/ A.U.)/ (km/CUL) 

CUL/astronomical unit 

POOL 

TABLE(45) 

(sec/ day)/ (sec/CUT) 

CUT/day 

POOL 

TABLE(46) 

(sec/hr)/ (sec/ CUT) 

CUT/hr 

POOL 

TABLE(47) 

(km/ CUL)/ (sec/ CUT) 

Convert CUL/CUT to 
kilometers/ second 

POOL 

TABLE(48) 

(km/CUL)/(sec/CUT) * (sec/hr) 

Convert CUL/CUT to 
kilometers/hour 

POOL 

TABLE (49) 

in rad/ sec * (sec/CUT) 

earth rotation in rad/CUT 

POOL 

TABLE(50) 

(km/ CUL)/ (km/ mi) 

miles/CUL 

POOL 

TABLE(51) 

(mi/CUL) * (CUT/hr) 

Convert CUL/CUT to miles/hour 

POOL 

TABLE(52) 

2 f - f 2 

e 2 - eccentricity of the earth 
squared 

POOL 

TABLE (53) 

e 

e - eccentricity of the earth 

POOL 

TABLE(54) 

0.0 

X component of U2 

POOL 

TABLE(55) 

COS € 

Y component of U2 

POOL 


155 





/BPOOL/ (continued) 


Variable 

Value 

Description 

Program 
Where Defined 

TABLE(56) 

sin e 

Z component of U2 

POOL 

TABLE(57) 

rate of change of the mean long, of sun in 
rad/C UT 

f - rate of change of the mean 
longitude of sun in rad/CUT 

POOL 

TABLE(58) 

J 2 * R 2 * 3/2 

J 

POOL 

TABLE(59) 

(sec/ CUT>/60 

minutes/ CUT 

POOL 

TABLE(60) 

(km/CUL)/(sec/CUT) * 1000 

Convert CUL/CUT to 
meters/second 

POOL 

TABLE(61) 

J 5 * R 2 /2 

Kj 

POOL 

TABLE(62) 

-Jj * RJ 

x 3 

POOL 

TABLE(63) 

-J 4 * R 4 * (3/8) 

k 4 

POOL 

TABLE(64) 

-J 5 * R 5 

K s 

POOL 

TABLE(65) 

(deg/hr)/ (deg/ rad) 

radian/hour 

POOL 

TABLE(66) 

(meters/ft) * 1000 

kilometers/ foot 

POOL 

TABLE(67) 

.40915751 

e - obliquity of ecliptic in radians 

POOL 

TABLE(68) 

1.0 

KMULT for drag 

POOL 

TABLE(69) 

in rad/ sec * (deg/rad) * 60 

earth rotation in degrees/minutes 

POOL 

TABLE(70) 

l/(min/CUT) 

CUT/ minute 

POOL 

TABLE(71) 

1.0 

X component of U1 

POOL 

TABLE(72) 

0.0 

Y component of U1 

POOL 

TABLE(73) 

0.0 

Z component of U1 

POOL 

TABLE (74) 

0.0 

tolerance for R • U 

POOL 

TABLE(75) 

1.0 

tolerance for magnitude RXU 

POOL 

TABLE(76)-(80) 

0.0 

not used 

POOL 


/DAFPRM/ 

COMMON/DAFPRM/LDAF(13) 


Variable 

Description 

Program 
Where Defined 

LDAF(13) 

Intermediate values from BRWORB used 
in computing the Data Acquisition Facility 
Parameters 

BRWORB 
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LDAFi = 


LDAF 2 


ldaf 3 


ldaf* 

4 


LDAF S = 


LDAF 6 = 


_1 

8 


y\ e" v 2 [1 " H£ 2 - 40£ 4 (1 - 50 2 )- 1 ] 


_5_ 74 
12 


e" r, 2 [1- 30 2 - 8<9 4 ( 1 - 50 2 )' 1 ] 

?2 



73 2 • 

~ 77 sin I 

?2 


ri 



64 ~ 2 Sin r ’ (4 + 3e " 2 > [1 " 902 " 24£ 4 (! “ 5 6 2 )- 1 ] 

7 2 V 


35 7s 

384 -y 1 
' 2 


e" 2 77 2 sin I" [l - 50 2 - 160 4 (1 - 50 2 )- 1 ] 


'V * 


y 2 r 


[1 - 116> 2 - 40£ 4 (1 - 5<9 2 ) -1 ] 


l|T7 ^ [1 - 30 2 - 80 4 (1 - 5 <9 2 )~ 1 ] 

/ 7 



Zi 

7; 


sin I 


• I 


_5_ zkul 
64 y' 2 e" 


sin I" (4 + 9 e" 2 ) [1 - 9^2 _ 2 4 i 9 4 ^ j _ 552^-1] j 


35 7 5 


384 


7, 


77 3 e" sin I" [1 - 50 2 - 160 4 (1 - 5<9 2 )-l] 
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ldaf 7 = 


LDAF 8 = 


LDAF 9 = 


' 16 y 2 [ + (2 + e" 2 >- 11 ( 2 + 3e" 2 ) 9 2 - 40 (2 + 5e" 2 ) <9 4 ( 1 - S0 2 )~ ] 


4OOe" 2 0 6 (1 - 50 2 )- 2 ] +4: — [2 + 

'V 1 


e" 2 -3(2 + 3e" 2 'i 0 2 


-8(1+ 5e" 2 ) 0 A (1 - 50 2 )- 1 - 80e" 2 9 6 (1 - 50 2 )" 2 ] 


J_1 _ / sin I" e" (92 \ 5 7s 

l 4 y ' 2 \ e" sin I"/ + 64 


/ 7] 2 sin I 
\ e" 


e" 9 2 
sin I" 


(4 + 3e" 2 ) + e" sin I" (26 + 9e" 2 ) 


[1 - 90 2 


- 24 0 A ( 1 - 50 2 )- 1 ] 


15 ^5 

32 


e" 6 2 sin I" (4 + 3e" 2 ) 


[3 + 169 2 ( 1 - 5 <9 2 ) _ 1 + 40£ 4 ( 1 - 50 2 )- 2 ] 


35 y s 


1152 


7 


2 L 


e sin 


I" (3 + 2 e" 2 ) - — 


" 3 e 2 


sin I' 


[1 - 50 2 


35 


166,4 (1 - 55 2 )- 1 ] + •=■ — e " 3 6 2 sin I" [5 + 326> 2 (1 - 5 6 2 )~ l 


y 


+ 80 <9 4 (1 - 56> 2 )- 2 ]j 
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LDAFjo = |- \y' 2 e '' 2e t 11 + 80# 2 C 1 - 5# 2 )" 1 + 200 # 4 C 1 ~ S# 2 )" 2 ] 


+ e " 2 9 [3 + 1602 (1 - 50 2 )- 1 + 40# 4 (1 - 5# 2 )- 2 ]} 

12 7 2 J 


LDAF n 


J_L zi e " g + j. _zi e " 6 

\ 4 y 2 s ^ n I" 72 s ^ n I 


(.4 + 3e" 2 ) [1 - 90 2 


- 24 # 4 (1 - 5 <9 2 )~ 1 ] + -|| —7 e" e sin I" (4 + 3e" 2 ) 

72 


[3 + 16# 2 (1 - 5 # 2 )- 1 + 40# 4 (1 - 5# 


2 >- 2 ] 


ldaf 12 


35 e" 3 6 

1152 sin I' 


[1 - 5<9 2 - 16# 4 (1 - 5# 2 )- 1 ] 


35 

576 v 

' 2 


e" 3 9 sin I" [5 + 32# 2 (1 - 5# 2 )" 1 


+ 80# 4 (1 - 


5# 2 )- 2 ] 


ldaf 13 


r ) 2 tan I 
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/DELKEP/ 

C0MM0N/DELKEP/DKEP(6) 


Variable 


DKEP(6) 


Description 

Perturbation A a, Ae, A i, A g, A h, A 1 


AX 


x t0 - X. 


'obs 


Program 
Where Defined 

BRWORB 


where 

X t0 = value of X on the pert tape at epoch 

X obs = value of X on the pert tape at the 
requested time 


Variable 

NPT 


/PIND/ . 

C OMMON/ PIND/N PT 


Description 


Perturbation tape logical number 

> 0 read pert tape 
< 0 do not read pert tape 
= 0 error 


Program 
Where Defined 

MAIN - Initializes 

ELEMLD 
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/PRTKEP/ 

COMMON/ PR TKEP/PKEP(3) 


Variable 

Description 

Program 
Where Defined 

PKEP(l) 

g 0 + Ag 

BRWORB 

PKEP(2) 

h 0 + Ah 


PKEP(3) 

1 0 + A1 



where 



X 0 = value of X at epoch 
AX = X t0 - X obs 



X t0 = value of X on the pert tape at epoch 


X obs = value of X on the pert tape at the 
requested time 


/RADIAN/ 



COMMON/RAD IAN/ TAU, AMBDA 


Variable 

Description 

Program 
Where Defined 

TAU 

t - longitude {in radians) of the sun on 
reference date 

MAIN 

AMBDA 

X - hour angle (in radians) of the first 
point of Aries on reference date 

MAIN 


2. BULLETIN COMMON Block Cross Reference Table 

This section contains a cross reference table describing the COMMON area 
structure in the BULLETIN program. 
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COMMON BLOCKS 


BULLETIN COMMON CROSS REFERENCE TABLE 


SUBROUTINES 


DDEEGHI J JKLN 


R S S S T V W 

R H C A DM RRLLTMNDUEAO PE PPEDPS IUUEM 

MWAK tdasdaecerstslpgdnarpprrdfapmtvcaz 

AONS DAAYTRGFOMATPCHLRAS GTOQE I UWCTEG I TPE 

I rpuan tcoalonlcolumrnlp efoudncoeh ts jolr 

N B L M F 0 ETRGDD0 DER0 TS 10 XT I0LVS TEELTBTKRDO 


NBLMF0 ETRGDD0 DER 


BPOOL XXX 

DAFPRM X 

DELKEP X^™ 

DHMNOD _X 

DHMNS 

DOTELE _X_ 

ETAP X 

GMPR X ; 

LPPRM X__ 

NOD X 

NODMAP _X 

NSMAP 

OSCELE X 

PERTL X 

PIND _X__X 

PRTKEP X 

RADIAN _X 

SOT 

SECPRM ~~X 

THETAP X 


X X X X 


X X X X X X 


X X X X X X 


>-i>o ^norni 


D. ABBREVIATIONS 


The following abbreviations and symbols are used in this program. 


a 

B 

CUL 


CUT 

e 

E 

f 


g 

h 

i 

K 

fi 


AC 

M 

N(2,Q) 

N(3,Q) 

P 

a 

p 

a 

r 

T o 

T(P,Q) 

X ^ 

Y > 
z 


Semi-major Axis 

Polar Radius of Earth 

Canonical Unit of Length 
(CUL = earth radius) 

Canonical Unit of Time 

Eccentricity 

Eccentricity Anomaly 

Flattening Coefficient (1/295.25) 

Brouwer's Notation for Argument of Perigee 

Brouwer's Notation for Right Ascension of Ascending Node 

Inclination 

Brouwer's Notation for Zonal harmonics 
Brouwer's Notation for Mean Anomaly 
Correction to Brouwer's Mean Anomaly 
Mean Anomaly 

First Order Drag Coefficients 
Second Order Drag Coefficients 
Anomalistic Period 
Derivative of Anomalistic Period 
Satellite Position Vector 
Unit Satellite Position Vector 
Epoch — Time of Elements 
Time of Drags 

Satellite Position Vectors 
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Satellite Velocity Vectors 


Right Ascension of Ascending Node 

Derivative of Right Ascension of Ascending Node 

Argument of Perigee 

Derivative of Argument of Perigee 

Rotation of Earth 

True Anomaly 

Hour Angle (in radians) 

Longitude (in radians) of Sun on Reference Date 
Motion of r 
Elevation 
180° in radians 

Gravitational Constant x Mass of Earth (Usually = 1) 



E. REFERENCES 


1. Method of Orbit Determination — Escobal, Pedro Ramon: John Wiley and 
Sons, Inc., 1965. 

2. Mathematics of Orbit Determination — Siry, Joseph W.: GSFC Publication 
X-547-64-151, June 1964. 

3. Brouwer- Lyddane Orbit Generator Routine — Galbreath, E. A.: GSFC 
Publication X-553-70-223, June 1970. 

4. Definitive Orbit Determination System: Module Performance and Design 
Descriptions, Volume II — IBM Corporation, under NASA Contract No. 

NAS 5-10022, November 1968. 
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BROUWER BULLETIN ROUTINE 


III. OPERATING INSTRUCTIONS 


This routine provides an economical means of disseminating pertinent spacecraft 
orbital information to observing stations and other interested persons. The 
Bulletin information includes the general characteristics of the orbit of the 
satellite, revolution numbers, as well as data useful for certain prediction pur- 
poses. An ephemeris is furnished to those who utilize the spacecraft for scien- 
tific, technological, and other purposes. Sufficient information from which the 
pointing angles may be determined from standard transformations is also 
included. 

The Bulletin routine reads input data from cards and, optionally, from the com- 
plementary perturbation tape. An output tape is produced containing pertinent 
spacecraft orbital information. 


A. REQUIREMENTS AND OPTIONS 

A production run requires a program tape containing the Bulletin and ancillary 
routines in the object module. The program tape is compiled using the appro- 
priate source deck and the Fortran H compiler of an IBM S/360 model 75 or 
model 95. Other requirements include three 9-track tape drives or two 9-track 
tape drives and one 7-track tape drive, a card reader, and an on-line printer. 
Data cards are prepared as specified in Section III B-2 (Input Card Format). 

Options are available in the program to load drag data from cards, to use the 
complementary perturbation tape, and to add the data acquisition facility 
parameters. In addition, this program has a change of constants capability. 

By use of the change of constant cards, any of the eighty values in TABLE can 
be changed if the user desires. The first forty values in TABLE must be 
changed by the first change of constant card and the last forty values by the 
second change of constant card. TABLE(80) is listed in COMMON/BPOOL/ in 
Section II C-l. 

There are special options which may be indicated by the change of constant 
cards. 

1. The sator output is suppressed by setting TABLE(36) equal to +1. 

2. The radio frequencies are stored in TABLE(37), (38) and (39). These 
locations are now set to zero but any or all may be changed. 
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3. The SPACEL osculating epoch elements output is suppressed by setting 
TABLE (40) equal zero. 

4. When creating the perturbation tape, drag can be included in the com- 
putation of the elements . Drag data can also be loaded from card to 
provide corrections to the Brouwer elements. The Bulletin user may 
desire to use drag from card only if no drag is on the perturbation 
tape or he may desire to use drag from tape and card. This option is 
controlled by TABLE(68), delta M drag multiplier. TABLE(68) is 
presently set equal to +1, which causes the drag data from card to be 
used whether or not drag is included on the pert tape. If TABLE(68) 
is set to 0, drag data from card is used only when drag is not included 
on the pert tape. If KMULT, which is on the pert tape equals +1, no 
drag is included on tape, and if KMULT equals 0, drag is included. 


B. INPUT 

Two types of input data are provided. 

1. Fixed — formatted cards are used to input values of the epoch elements, 
observation times and various options. 

2. The complementary perturbation tape provides corrections used by the 
Brouwer Orbit Generator; its optional use is controlled by the PERT 
option on the elements time card. 

1. Limitations 

There are two limitations on the card input. 

a. 27 is the maximum number of columns to be computed from a single 
column card (see option 1 for card h). If the user inputs a number 
larger than 27, an error will not be generated but the number will be 
reduced to 27. 

b. For SPACEL elements output, column elements time must be equal to 
a multiple of 0.05 day. 

2. Input Card Order 


a. Run identification 
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b. Change of constant card(s) followed by blank, or blank if no constant is 
to be changed 

c. Satellite identification 

d. Drag card(s) followed by blank, or blank if no drag is used 

e. Change of constant card(s) followed by blank, or blank if no constant is 
to be changed 

f. Elements time (Epoch) 

g. Elements (2) 

h. Column elements time and drag data card(s) followed by a blank. 

Option 1 — 1 or 2 cards specifying the number of columns and the 
delta T desired for each card 

Option 2 — 3K + 1 cards spaced at equal time interval, K = 1 to 18 

i. Date of issue 

j. Remarks 

k. Bulletin request 

l. Spacel bulletin 

n. Data Acquisition Facility (DAF) parameters card or blank if not re- 
questing DAF parameters 

3. Card Format 


Format Specification Interpretation 

Key: n,m = integer numbers 

b = a blank space 

Format Code Interpretation 

In digits with no decimal point right adjusted in a field of 

n columns 

example: 35 in an 13 format: b35 
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Format Code 
Fn. m 


Dn. m 


An 


Interpretation 

digits with a decimal point anywhere in a field of n 
columns or digits punched with no decimal point, in 
which case the point will be assumed between the 
m - 1th and mth column of the field 

example: 40.1 in an F5.2 format: b40.1 or 40.1b 
or b401b 

digits with a decimal point anywhere in a field of n 
columns or digits right justified in a field of n columns 
with an exponent of the form D ± XX, where XX is the 
power of 10 to which the number is raised. If there 
is no decimal point it will be assumed to be m places 
to the left of the "D". 

example: 40.1 in a D8.5 format: b40.1bbb or 
40.1 D + 00 or 401. D - 01, etc. 

n alphanumeric characters 


a. Run Identification Card 



00000000 
1 7 3 4 9 8 7 8 
II 1 1 11 II 


0 0 0 0 0 0 0 0 0 0 0 0 0 0 g G 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

8 18 II 17 13 <4 13 IB 17 18 19 30 31 33 33 34 33 38 33 38 38 30 3> 33 33 34 33 38 37 38 39 4^ 

1 I I I I 1 I 1 I 1 I I I I I I I I 1 I 1 It 


72227222 
3 3 3 3 3 3 3 3|3 
44444444 

55555555 


ooooooooooooooooooooo ooooooooolo 

41 43 43 44 43 48 4 7 48 48 30 31 33 33 34 33 38 37 31 39 80 81 S3 83 84 83 88 87 88 89 70 

ill II I It I | | | | | | | | | | | 1 | 


RUN IDENTIFICATION CARD 


2222 2 272222222222222 2 2 2‘rTTTTTT 
3333333333333333333333333333333 
44444444444444440444444444444444 


000000000 
7113 73 74 79 7* 77 78 79 80 

I II I 1 I 1 11 I 



4 4 4 . 4 4 . rrt 2 2 2 2 2 2 2 2 7 7 2 2 2 2 2 7 2 2 2 7 

3333333333333333333333333333333 

0444444444444444444444444444444 

5 5 5 5 5 5 5 — ^«<i5 5 5 5 5 5 5 5 5 5 5|5 


7222727727 

333333333 

444444440 

5! 


Card Col. 

Format 

Field Description 

1-8 

bbb . . . b 

Leave blank 

9-40 

A40 

First line of run identification 

41-70 

A40 

Second line of run identification 

71-80 


Not used 
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b and e. Change of Constant Card 



|0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Q C 0 

3 4 1 s 7 a 9 10 n 17 13 l« 15 II 17 It n 2a 71 22 23 24 25 1 

I | 11 11 U 1 11 1 11 11 t H 1 1 1 1 


0 0 
« 2 
1 1 

2 2|2 222222222222222222222 2[2 2 

3 

4 4 

5 


3|3 333333333333333333333 3 
44444444444444444444444 
515 5 5 5 5 S 


0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
III M 30 3< 33 33 34 35 X 37 31 31 40 4) « 13 44 4S4I47 (3 43 30 3 


CHANCE OF CONSTANT CARD 


1 1 


2 2 


3333333333333333333333333 


| 4 4444444444444444444444 
5555555555555 


2 2 


44 


0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
|>3 m u m o) u si oo (i 12 u u u k u tt n> ?i n n m ni 

1 I I I I I 1 II I I I 1 I 11 1 M I 1 t if 
2222222222222222222222 


333333333333333333333333 


4444444444444444444444 4 | 

5555555555555 


Card Col. 

Format 

Field Description 

1-2 

12 

Table location to be changed 

3-25 

D23.16 

Constant to be inserted in above location 

26-27 

12 

Table location to be changed 

28-50 

D23.16 

Constant to be inserted in above location 

51-52 

12 

Table location to be changed 

58-75 

D23.16 

Constant to be inserted in above location 

76-80 


Not used 
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c. Satellite Identification Card 



0 0 0 0 0 0 0 0 0 0 0 0|0 0 0 0 0 0 0 0 0 0 0 0|0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

|si si si w (i uumukoiuu » n nn )< »i» n hum 

1 1 1 1 1 1 1 I 1 t 1 1 I 1 I I I I I 1 1 I 1 I 


IpJ 34 3S 01 11 30 il 40 41 42 42 44)45 45 41 44 45 50 51 52 53 54 55 54 

inn 


SATELLITE IDENTIFICATION CARD 


1333333333333 
1444444444444 
1555555555 


33333333333 3 | 
|4 44444444444 
5 5 5 5 . 


2 2 2 2 2|2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 
333333333333333333333333 

44444444444441444404444O 

5555555555555! 


Card Col. 

Format 

Field Description 

1-7 

17 

Satellite identification number 

8-13 

16 

Date of reference: year, month, day 

14-22 

19 

Hour angle of the first point of Aries on 
above date (k): hours, minutes, seconds 

23-32 

no 

Longitude of Sun (r): degrees, minutes, 
seconds 

33-44 

b 

Leave blank 

45-56 

A12 

Satellite name 

57-80 


Not used 
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d. Drag Card 



2 2 2 2 2 2 2 

3 3 3 3 3 3 3 

4 4 4 4 4 4 4 
5555555 


2 2 


2 2 


3 3 3 3 


0 0 0 0 0 3 0 0 0 0 0 0 0 0 0 0 0 0 B 0 0 0 a 

i|24 2; 22 29 2D Ji 22 22 24 22 21 22 24 2) 40 *1 42 42 44 42 44 41 44(4! 

11 1 I I 11 I I .1 I 1 I 1 

DRAG CARD 

2 2 2 2 2 2 2 2 2 » 2 2 2 2 ! 


3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 Bj 
44444444444444444444444 
5555555555555 



OODOOODOOOOOOOOOOOOOOOO 

12 50 21 52 23 24 22 24 22 5) 5)40 42(2 (3 44 42)4 02 4)4) 30 2I| 
11111111111111111111111 


22222222222222222222222 
33333333333333333333333 
44444444444444444444444 
555555555555 


000000000 
22 33 24 12 14 32 34 3)10 
111111111 

2 2 2 2 2 2 2 2 2 

333333333 

|4 44444444 

5J 


Card Col. 

Format 

Field Description 

1-7 

17 

Satellite Identification number 

8-13 

16 

Date of Drag: year, month, day 

14 

b 

Leave blank 

15-18 

14 

Time of Drag: hours , minutes 

19 

b 

Leave blank 

20-24 

15 

Seconds 

25 

b 

Leave blank 

26-48 

D23.16 

N(2,Q) First Brouwer drag parameter 

49-71 

72-80 

D23.16 

N(3,Q) Second Brouwer drag parameter 
Not used 
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f. Element Time Card 



2 2 2 2 2 2 2 
3 3 3 3 3 3 3 
4444144 
5 5 5 5 5 5 5 


2 2 2 
3 3 3 

14 4 4 


OOOOOOOOOOOOOOOOOOQOOOOOOOOOOOOOOOOOO 

||K II IS IS 10 ]t 3; 33 33 33 3S 33 30 31 3D 41 33 33 33 33 SS 3) 3| 31 50 SI S3 53 S3 S3 S3 31 SI S3 SO SI S3 

iniii ; 1111111I1111M1 


2 2 2 2 2 2 


ELEMENT TIME CARD 


>222222222222222 


3333333333333333333333333333333333333 
44444444444444404444444444444444444441 

J 5 5 5| 




00000000000 
lonnnuisisnisisso 
1 1 I I 1 1 I 1 1 I I 

2222222222 ? 

33333333333 

44444444440 

55! 


Card Col. 

Format 

Field Description 

1-7 

17 

Satellite Identification Number 

8 

b 

Leave blank 

9-14 

16 

Date of elements: year, month, day 

15-18 

14 

Time of elements: hours, minutes 

19-24 

F6.3 

Seconds 

25 

11 

Input element type code: 

1 = type 1 elements on next 2 cards 

2 = type 2 elements on next 2 cards 

26-62 

bbb . . . b 

Leave blank 

63-67 

15 

Pass number of elements 

68 

b 

Leave blank 

69 

11 

Perturbation tape option 
1 = pert tape is used 
blank = no perturbation 

70-80 

Not used 
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g. Elements Card 1 



0 0 D 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

1 2 HU 2 1 9 to It 12 13 14 K tt 1? HU 20 21 22 23 
11111111111111111111111 

22222222222222222222222 

33333333333333333333333 

44444444444444444444444 

55555555 


0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 D 0 D 0 

|»24 21 30 23 30 31 32 33 34 33 30 33 30 33 40 41 42 43 44 43 44 41 

11111111 i ■ — - j 1 

22222222 


ELEMENTS CARD 1 


333333333333333333333333 


04444444444444444444444 

5 5 5 5 5 5 5 5 5 5 5 5 


JOOOOOOOOOOOOO 0 0 0 0 0 0 0 0 0 0 

|4fl|lS 30 31 31 33 34 33 50 31 31 5| SO SI II 13 34 S3 04 01 40 60 10 It 

1 1 1 1 1 1 1 1 1 1 1 1 1 I 1 I I 1 1 1 I 1 I 
22222222222222222222222 
33333333333333333333333 

144444444444444444444444 
555555555555 



000000000 
I3i3i4!3nn is nao 
111111111 

22222222 ? 

333333333 

444444440 

5 J 


Card Col. 

Format 

Field Description 

Type 1 

Type 2 

1-23 

D23.16 

X (CUL or km) 

a (CUL or km) 

24 

b 



25-47 

D23.16 

Y (CUL or km) 

e 

48 

b 



49-71 

D23.16 

Z (CUL or km) 

i (radians or degrees) 

72-80 


Not used 

1 
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g. Element Card 2 



0 0 0 0 0 0 0 0 D 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

1 2 3 * 1 * 2 ■ 9 10 II 22 13 I* IS K U ll 19 20 21 22 23p4j 

t I 11 1 1 1 I 1 1 1 1 I 1 I 1 I 1 I t 1 1 l| 1 


22222222222222222222222 
33333333333333333333333 
44444444444444444444444 
55555555 


0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

lilt 11 29 29 10 91 32 31 34 39 39 32 30 39 40 41 42 43 4* 95 II «j| 

I I t | 1 I I I , ,• I 


ELEMENTS CARD 2 

2122222222 hrrTTTrn rrri 2| 


33333333333333333333333 3! 


44444444444444444444444 


5|B 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

[*S 20 91 92 93 94 33 99 9) 91 99 40 91 92 99 94 43 94 02 99 41 20 2|| 

1 1 1 1 1 I I 1 I I 1 I 1 1 1 1 1 I I 1 I 1 


2222222222222222222222 
3333333333333333333333 
| 4 4444444444444444444444 
555555555555 



000000000 
72 93 2* 29 29 21 24 34 90 

1 1 1 n n n 


2 2 2 2 2 2 2 2 2 2 
3333333333 

444444440 

5 ■ 


Card Col. 

Format 

Field Description 

Type 1 

Type 2 

1-23 

D23.16 

X (CUL/CUT or km/ sec) 

M (radians or degrees) 

24 

b 



25-47 

D23.16 

Y (CUL/CUT or km/ sec) 

u) (radians or degrees) 

48 

b 



49-71 

D23.16 

Z (CUL/CUT or km/ sec) 

Q (radians or degrees) 

72-80 


Not used 

1 


Cards 1 and 2 will have one of the above two formats (type 1 or type 2) depending 
on the type specified in column 25 of the elements time card. 

Elements on card 1 and card 2 must be of the same units (CUL - radians or 
km - degrees). 
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Column Elements Time and Drag Data Card 


COLUMN N(2,Q) 


COLUMN N (3, Q) 


0000000 

1 2 3 * 5 6 1 
1111111 

2 2 2 2 2 2 2 

3 3 3 3 3 3 3 

4 4 4 4 4 4 4 

5 5 5 0 5 5 5 


0 0 0 0 0 0 0 C 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

1) II It to 31 31 33 2* 23 26 21 21 39 30 31 22 33 34 33 31 31 31 31 40 41 42 43 44 43 41 4 ) 41 43 30 31 S? 33 54 33 36 31 31 33 60 61 62 63 34 63 66 03 61 61 

1 1 1 II 1 1 1 1 1 1 I 1 1 ' ' * .- 1 .. '- | 1 1 1 I 1 1 1 1 1 I 1 1 1 1 1 1 1 1 I 1 I 

COLUMN ELEMENTS CARD 

2 2 2 2 2 2 2 2 2 2 2 2 2 2 177777777777777771 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 

3 3 f 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 

44444444044444444444444404444444444444444444444444444 

"" “ ,,, **Ai5 5 5 5 5 5 5 5 5 5 5 5 5>^ _ 5 5 5 5 5 5 5 5 5 5 



Card Col. 

Format 

Field Description 

1-7 

17 

Satellite Identification Number 

8 

b 

Leave blank 

9-14 

16 

Date of column elements: year, month, day 

15-23 

19 

Time of column elements: hours, minutes, 
seconds 

24-46 

D23.16 

N (2,Q) First column drag parameter 

47-69 

D23.16 

N (3 ,Q) Second columns drag parameter 

70-71* 

12 

Number of columns to be computed 

72-78 

17 

Column delta T in minutes 

79 

b 

Leave blank 

80 

n 

Punch if two cards are to be read. 
Blank if only one card is to be read. 


Note: For spacel elements, column time must be equal to a multiple of 0.05 day. 

* Columns 70-80 are optional. 

Option 1 — If columns 70-80 are used, one or two cards are needed, specifying the number 
of columns and delta T desired for each card. If two cards are used, column 80 
of the first card must be punched. 

Option 2 — If columns 70-80 are not used, 3K + I column cards, spaced at equal time inter- 
val, must be loaded (K = 1, 18). 
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Issue Date of Bulletin Card 



0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Q 0 Q 0 0 0 0 0 C 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Q 0 0 0 D U D 0 U 0 U 0 0 U U U U U u u u u u u u u u u u u u u 

1 2 1 4 * S J I 3 10 11 13 >3 14 'MS H II i3 20 21 22 33 34 35 26 31 31 29 30 Ji II 33 31 35 31 33 31 39 40 II <2 43 44 IS «6 41 46 41 50 51 53 53 54 96 M 51 51 5160 61 12 63 64 IS M 61 61 11 10 11 13 73 74 15 15 17 II It 10 

1| M M 1 1 1 I M U 11 I 1 1 11 1 11 11 1 1 11 t I 1, 1 I 11 I I I I I 1 I 1 II 1 1 II 1 1 1 I 1 1 • I 1 ' 1 11 I I > 

ISSUE DATE CARD 

2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 W-rrrr^T^-r 4 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 

3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 

4444444444444444444444444444444444444444444444444444444444444444444M0444044**** 


5 5 5 5 5 5 5 5j 


l5 5 5 5 5 5 5 5 5 5 5 5 5J 


.5555555555555! 


Card Col. 


Format 


Field Description 


Issue Date of Bulletin: year, month, day 
Not used 


j. Bulletin Remarks Card 


LINE 1 OF REMARKS 



LINE 2 OF REMARKS 


0 0 0 0 0 0 0 0 0 0 0 0 D 0 0 0 0 0 D 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

35 36 31 36 31 40141 43 43 44 45 46 41 46 45 50 51 5? 53 94 15 56 II 51 51 60 61 63 63 64 65 16 01 66 66 10 n 13 

’1 ’ ' ’ ’ 1 I I I I II II I I I I II I 11 1 I 1 1 11 M I 

REMARKS CARD 

.......... . V 222222222222222222222222222 


33333333333333333333333333333333333333 




Card Col. 

Format 

Field Description 

1-8 

bbb . . . b 

Leave blank 

9-40 

A40 

Line 1 of Remarks 

41-72 

A40 

Line 2 of Remarks 

73-80 


Not used 


177 














k. Bulletin Request Card 


WMAP 

1 

START 

DATE 

0 0 0 0 

0 0 

0 0 0 0 0 0 

17 10 

S 6 

10 9 10 II 17 

llll 

1 1 

11 11 1 I 

2 2 2 2 

2 2 

2 2 2 2 2 2 

3 3 3 3 

3 3 

3 3 3 3 3 3 

4 4 4 4 

4 4 

4 4 4 4 4 4 

5 5 5 5 

5 5 

5_5j>fc rf, ‘ 



START 

TIME 


END 

DATE 


END 

TIME 

0 

0 0 0 0 0 0 

0 

0 0 0 0 0 0 

0 

0 0 0 0 0 0 C 

13 

1 

U IS IB I? ta 13 

1 1 11 1 1 

20 

1 

21 2? 23 24 23 ft 

111111 

21 

2S 29 30 3i 22 33 3 

111111 

2 

2 2 2 2 2 2 

2 

2 2 2 2 2 2 

2 

222222. 

3 

3 3 3 3 3 3 

3 

3 3 3 3 3 3 

3 

3 3 3 3 3 3 3 

4 

4 4 4 4 4 4 

4 

4 4 4 4 4 4 

4 

4 4 4 4 4 4 4 


INC. 


I|0 0 0 0 0 0 0 0 0 0 0 0 0 0 G 0 0 0 0 0 0 0 D 0 0 0 0 0 0 0 0 C 3 0 0 0 0 0 0 0 0 0 0 

ijis 40 41 40 4]|*4 45 4S 4J 40 4S SO SI 57 S3 54 SS SO SI SI SS 00 61 01 (3 04 IS K 01 U 41 10 Tl 77 IS 14 7S 10 n 10 II 10 

1 I 1 1 1 1 II I I I I I 1 1 It 1 1 1 1 1 1 I I I I I II II I I I 


REQUEST CARD 



22222222222222222222222222222222222 
3333233333333333333333333333333333333 

444444444444444444444444444444044044O 
5 5 5 5 5>*^ 5 5 5 5 5 5 5 5 5 5 5 5 5 ! 


Card Col. 

Format 

Field Description 

1-4 

A4 

The word WMAP 

5-6 

bb 

Leave blank 

7-12 

16 

Starting Date: year, month, day 

13 

b 

Leave blank 

14-19 

16 

Starting Time: hours, minutes, seconds 

20 

b 

Leave blank 

21-26 

16 

Ending Date: year, month, day 

27 

b 

Leave blank 

28-33 

16 

Ending Time: hours, minutes, seconds 

34 

b 

Leave blank 

35-37 

13 

Latitude increment (degrees) 

38 

b 

Leave blank 

39-43 

15 

Inclination (.01 degrees) 

44 

b 

Leave blank 

45-49 

50-80 

15 

Pass number 
Not used 
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1. Spacel Bulletin 



Card Col. 

Format 

t — — 

Field Description 

1-7 

17 

Satellite Identification Number 

8 

b 

Leave blank 

9-14 

A6 

Abbreviation of popular name of spacecraft 

15 

b 

Leave blank 

16 

A1 

Source of spacecraft 

17 

A1 

Launching site 

18 

b 

Leave blank 

19-22 

A4 

Norad number of object 

23-29 

bbb . . . b 

Leave blank 

30-37 

F8.3 

Initial Julian Date of Space associated with 
the first time the spacecraft achieved its 
orbit 

38-44 

b 

Leave blank 

45-48 

A4 

Brightness 

49 

b 

Leave blank 

50-54 

A5 

Ratio of spacecraft weight to reference 
weight in kilograms per square meters 

55 

b 

Leave blank 

56-60 

A5 

Radio frequency in deeakilocycles 

61 

A1 

Type of modulation 

62-64 

A3 

Transmitted Power in Centiwatts 

65-66 

b 

Leave blank 

67 

A1 

Center of Attraction 

68-69 

A3 

Orbit number 

70-80 


Not used 
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m. Data Acquisition Facility Parameter Request Card 




1 

TIME 

ll 

□ 

a 

0 

0 0 

0 0 

15 

t 

IE 17 

1 1 

11 19 

I 1 

2 

22 

2 2 

3 

3 3 

3 3 

4 

4 4 

4 4 


|20 21 22 23 24 25 28 21 28 28 30 3' 32 33 34 3 5 38 3 7 38 5 8 40 41 42 43 44 45 48 4 7 48 48 50 51 52 53 54 55 5* 57 81 58 SO (I (2 83 84 (5 M 0/ M 8* 70 21 72 73 M 35 78 71 78 U 

1 1 l t I n i I I I i t t i i i i i i i, i tiiiiiiininiiiiMiMiniiiiiMiiM 


2222212222222222 


DAF CARD 


2222222222222222222222222222222222221 


444444444444444444444444444444444444444444444444444444a44444a 
5555555555555' 


Card Col. 

Format 

Field Description 

1-7 

17 

Satellite Identification Number 

8 

b 

Leave blank 

9-14 

16 

Date of observation: year, month, day 

15 

b 

Leave blank 

16-19 

14 

Time of observation: hours, minutes 

20-80 


Not used 


T M.JIJN — 5 BROUWER RULLETIN 
05+0. 80 681? <+7000000000+ 03 

6R06 60 271020108421 16633 112876600 IN JUM— 5 

6806602710220 0000-00000 +0 . 1 60 39 OOOOOOOOOOOD-OR 

6R06602 710220000000.0002 

* 0 . 12510 R4 511 94000 00 + 0 1 +0. 1 1 5 7 6 1 7002 2 30000D+00 +0. 14079 37 93054000 OD+O 1 
+0. 1 72 7337R 69 1R0 0000+01 + 0 . 606780704 1 520000000 1 +0. 34870792983300000+00 
6806602 710223000000000+0. 16039000000000000-08 

710303 

INJIIN-5 BROUWFR BULLETIN 

WMAP 710223 000000 710302 0R1500 010 04090 11292 
6806602 IMJIIN5 UW 3338 4914 3914 3977. R 4 ? 11262 


RUN ID 
CONST 
BLANK 
SAT ID 
DRAG 
BLANK 
CONST 
EPOCH 
ELEM 1 
ELEM 2 
031008000 
BLANK 
I SSUE 
REMARKS 
REQUEST 
SPACEL 
DAF 


Figure 2. Listing of Sample Input Data 
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4. Perturbation Observation Tape Format 

The complementary perturbation tape provides corrections used by the Brouwer 
Orbit Generator and the Bulletin Prediction. Its optional use is controlled by 
the PERT option on the elements time card (card f). A punch in column 69 in- 
dicates that the tape is to be used. The pert tape is a 9-track binary tape con- 
sisting of one header record followed by seven or more data records. 

See Section III-C for the job control language cards required in a production 
run of the Bulletin including perturbations. 


Table 1 

Format of Perturbation Observation Tape 

Header Record 


Word 

Number 

Word Contents 

0 

Fortran word count 

1 

Time increment — days 

2 

Month 

3 

Day 

4 

Year 

5 

Satellite Identification Number 

6 

Input semi-major axis — e.r. 

7 

Input eccentricity 

8 

Input inclination 

9 

Input right ascension of the ascending node — degrees 

10 

Input argument of perigee — degrees 

11 

Input mean anomaly — degrees 

12 

Input time from midnight — days 

13 

Input period — minutes 

14 

Number of records on tape excluding header and trailer records 

15 

Delta mean anomaly option indicator (KMULT) 

(1 — delta drag mean anomaly not computed on tape 
0 — delta drag mean anomaly computed on tape) 
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Table 1 (continued) 


Date Record 


Word 

Number 

Word Contents 

0 

Fortran word count 

1 

Time in seconds from epoch 

2 

a (semi -major axis) — e.r. 

3 

e (eccentricity) 

4 

i (inclination) - tt /2 < i < 77/2 

5 

AM (delta mean anomaly change from T 0 ) <f> < AM < 277 

6 

a> (argument of perigee) < co < 2n 

7 

n (right ascension of the ascending node) <p < £1 < 2 rt 


Sentinel Record 


Word 

Number 

Word Contents 

0 

Fortran word count 

1 

.9999999999999999 x 10 30 

2-8 

Irrelevant 


End of File after sentinel record. 


C. SET UP AND RUNNING PROCEDURE 
1. Requirements 

IBM S/360 model 75 or model 95, three tape drives including three 9 -track or 
two 9-track and one 7-track, on-line card reader, and on-line printer. 
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2. Tape Assignments 


Table 2 

Tape Assignments 


Tape Function 

Tape Description 

Read Only Tapes 

9T 


Program Tape* 


9T 

Unit 4 

Perturbation Tape (Optional) 

Generated Tape to be Saved 

7T or 9T 

Unit 9 

Bulletin Output Tape 


*The program tape contains the Bulletin and ancillary routines in the object module. Presently 
the tape contains two sequential data sets,*the source file and the object module which was 
compiled from the source coding using the Fortran H compiler of an IBM S/360 model 75 or 
model 95. 


3. Card Reader 

Reads input data cards and required JCL cards. 

4. On-Line Printer 

Indicates whether a constant has been changed or not and writes the initial con- 
ditions and error messages. (Refer to Section III D.) 

No special paper, loop, or board requirements. 

See Section in B for a detailed description of the data cards and perturbation 
tape required. 
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5. IBM S/360 Job Control Language Cards 


/ /G5DBMBUL JOB 


// EXEC LINKGO, REGION. G0 = 300K 


//LINK.SYSLIN 

// 

// 

// GO .FT04F001 
// 

// 

/ / GO . F T09 F00 1 
// 

7/ 

// GO .SYSUOOMP 
//GO.DATA5 00 


DO UNIT=2400-9,0ISP=(0LD,PASS) » LABEL= ( 2 ♦ BLP ) 
DCB= (RECFM=FB t LRECL=80»BLKSIZE=3200*DEN=2) , 
VOLUME=SER= 1 647 G 

OD UNI T=2400-9,DISP= (OLD, PASS ) ,LABEL=( , BLP) , 
DCR=(RECFM=VS»BLKSIZE=128»LP. ECL = 124»DEN = 2* 
VOL UME = SE R = XXXXX 

DD UN I T=2400-7 »DI SP = ( NEW »PASS ) t LABEL= ( ,BLP) , 
0CB=(RECFM=FA,BLKSIZE=120,TRTCH=ET,DEN=1 ), 
VOLUME=SER= 

DD S YSOUT = A 


/* 

// 


-- DATA CARDS ** 
** DATA CARDS ** 


D. OUTPUT 

The Bulletin output is generated on an on-line printer and a 9-track or 7-track 
tape which is assigned to tape unit 9. The initial conditions, changed constants, 
and error messages are printed on both the on-line printer and tape. 

1. On-Line Printer 

The on-line printer is used to write run identification data, initial conditions 
data, changed constants, and error messages. A sample on-line printout is 
shown in Figure 3. Tables 3 and 4 list normal statements and error statements, 
respectively. 


Table 3 

Normal Statements 


Statement 

Explanation 

No DAF Parameters Computed 

User does not desire to print DAF 

In This Run 

parameters and has blank input DAF card. 

TABLE (XX) = YYYY 

Specified location XX of TABLE now 


equals indicated value YYYY. 

End Bulletin Run 

Program has followed normal path. 
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TABLE! S» = 0.806812* 2O0CC0OQCO C 3 


ribs o.p. brouwer bllletin routine 

INJUN-5 BROUWER BULLETIN 


INJUN-5 BRCUWER EULLETIN 


10. NO. REF. DATE LAMBDA MMS TAU OMS SATELLITE 

6006502 71 2 1 a 42 1160 311 2S 26600 INJUN-5 

INPUT OUANTI TIES FROM CAROS 

E“OCH 71 2 20 0 0 0 


* 6 IS C MEG A C OMEGA M 

0.125I0845D 01 0.11S761700 00 0.140703790 01 0.17273379C 01 0.606780700 01 0.348707930 00 

CONVERTED QUANTITIES 

x V 2 XCOT YDOT ZDOT 

-0.501 191920 00 0.200557510 00 0.912065540 00 -O .8 460 338 9 C CO 0 .964604000— 0 1- 0 . 51 5028700 00 

DRAG EFFECTS T (P.Q) 6 (2,01 N 13. C) 

710220 0 0.160390 000-08 0.0 


COLUMN TABLE TIP. 01 M2, 01 NI3.Q) 

710223 C C. 160390 000-08 0.0 


M I » M IN* 10 0 
3 1008000 


F4RTM CONSTANTS ML POTATION RACIUS FLATNESS 

O.ICCCOCOOO Cl 0 .8883 36900—01 0.100000000 01 0.3 35209 190-02 

HARMONICS 

*2 13 K4 KS 

0.5412*0000-03 0. 256000000— 05 0.69000000 C—0 6 u .600000 OOD— 07 


J 1 

0. 162372C0D-C2 C.O 


0.0 


INPUT UNITS CUL ANO RADIANS 

* e i 

0. 1 26 1O04S1 1 9400000 01 0. 1 1 37* 1 7 0022 300 000 00 0 . 1 4079379305400 OOD 01 
S OMEGA C CMEGA M 

0.17273378691000000 01 0 . 6 C670O 7 04 1 52 00 00 0 01 0 .34870772903300000 00 


CONVERTEO QUANTITIES KM ANO CEGREES 

A E 

0.79796246971023010 04 0.11576170022300000 00 
S OMEGA C CMEGA 

0.98969169697134700 02 0.3476697343788582 0 03 


I 

0 .8066890 123632*240 
M 

O .19970402662174950 


02 

02 


BROUWER OAF PARAMETERS 
NO OAF PARAMETERS COMPLTED IN THIS RUN 

END BULLETIN PUN 


Figure 3. Sample On-Line Printout 
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Table 4 


Error Statements 


Statement 

Explanation 

Incorrect Elements SAT. ID. 

Input satellite ID does not agree with 
elements time satellite ID.* 

Wrong SAT. ID. on Drag Card 
Drag Load Unsuccessful 

Input satellite ID does not agree with drag 
satellite ID. 

Error in Column Elements 
Request Card 
Column Time Table Load 
Unsuccessful 

The format of the column elements time 
and drag data card is incorrect. When 
using option 2 for the format, columns 
70-80 must be used. 

Wrong SAT. ID. on Column 
Elements Card 
Column Time Table Load 
Unsuccessful 

Input satellite ID does not agree with 
column elements time satellite ID. 

Elements Type Equals Zero 

The element type on the elements time 
card must be either 1 or 2 to indicate the 
type of elements used. Refer to Sec- 
tion III B-3 for the format and description 
of the elements time card. 

Latitude Increment Exceeds 
Inclination 

The latitude increment can not exceed the 
inclination on the Bulletin request card. 

Tape Check on Perturbation 
Tape 

An uncorrectable error occurred when 
reading the perturbation tape. If tape is 
the correct tape to be used, rerun job. 

Perturbation SAT. ID. 

Input satellite ID does not agree with pert 
tape satellite ID. 

Incorrect Data to CHANPL 

The location of the constant to be changed 
is not in the range of the Table (0 thru 80) . 

Error in DAF SAT. ID. 

Input satellite ID does not agree with 
DAF satellite ID. 
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Table 4 (Continued) 


Statement 

Explanation 

Error, Year of Reference, 

Error occurred in subroutine DREFOD. 

XXXX, Greater Than 

The year of reference must always be 

Observation Year, XXXX 

less than or equal to the observation year. 

No SPACEL Bulletin Output 

The column elements time is not equal to 

For This Run 

SPACEL Elements Do Not Have 
Accuracy of Column Elements 

a multiple of 0.05 day. 


2. Output Tape Format 

The Bulletin tape output consists of the following sections: 

1. Title Page 

2. Initial Conditions 

3. Space Elements 

4. Spacel Bulletin 

5. Equator Crossings 

6. One Orbit Ephemeris 

7. Sator Code 

8. Data Acquisition Facility Parameters 

The Spacel Bulletin, Sator Code and Data Acquisition Facility Parameters are 
printed according to user option. The other sections are always printed. 

A listing of the sample BCD output tape is shown in Figure 4. The title page 
contains the satellite identification, start and end times for computing equator 
crossings and user remarks. 

Input values of some of the pertinent parameters appear in the initial conditions 
as listed in Figure 4a. If the PERT tape is used, the first line printed is 
"Complementary Perturbations". 
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Figure 4b contains space elements which indicate mean characteristics of the 
orbit at the epoch. The subsection entitled ’’Descriptive Space Elements" con- 
tains various derived quantities which are of interest to observers at the earth’s 
surface. "Prediction Space Elements" provides elements for use when approxi- 
mate satellite positions are needed. 

Osculating Space Elements and their Cartesian equivalents are listed in Figure 
4c. The Cartesian equivalents are obtained using the given values of GM and 
the value of the earth’s equatorial radius. The double prime elements at the 
epoch are also listed and the position and velocity vectors are computed in sev- 
eral units . 

Key space element information contained in Figures 4b and 4c is combined and 
furnished in condensed form in the Spacel Bulletin, Figure 4d. Note that the 
last digit of each line is a counter to be used by observation stations and other 
interested persons. The counter is the sum of all digits in its line modulo ten. 
The last line of the Spacel Bulletin contains osculating space elements at the 
first requested time of the prediction space elements. 

The subsection entitled "Equator Crossings" contains each ascending nodal 
crossing from the requested start time to the end time, its revolution number, 
its date and time, and its west longitude in degrees as indicated in Figure 4e. 

Figure 4f contains a one orbit ephemeris for the middle nodal crossing shown 
in Figure 4e. The ephemeris gives the satellite positions at regular intervals 
according to the requested latitude increment. Time is specified in terms of 
minutes from the time of the ascending nodal crossing. Longitude is given in 
degrees and decimal fractions. Height above the ellipsoid is given in terms of 
a decimal fraction of a kilometer. Times when the satellite is not in the earth’s 
shadow are indicated by means of an asterisk following the height. 

Figure 4g contains the Brouwer Data Acquisition Facility Parameters. 
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6806602 INJU 

FROM 710223 0 

TO 710302 81530 

[ N J UN-5 BROUWER BULLETIN 

Figure 4. Listing of Sample Output Tape 
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GODDARD SPACE FLIGHT CENTER 
GODDARD ORBIT BULLETIN 
MAR 3,1971 

SATELLITE INJUN-5 



SPACE 

ELEMENTS 


EPOCH 


YRMODYHRMMSS. SSSS 


CALENDAR DATE 


71 220 000. 0 

UT2W 

JULIAN DATE FOR SPACE 

4903. 00 

UT2M 

PERIOD, ANOMALISTIC 


118.289055 

MIN 

PERIOD DERIVATIVE 


-0.03807 

MICRODAYS/ DAY 

ECCENTRICITY 


C. 11669938 


INCLINATION 


80.66783 

DEG 

RIGHT ASCENSION OF 

ASCENDING NODE 

347.65948 

DEG 

RT.ASC.OF ASC. NODE 

DERIVATIVE 

-0.25138 

DEG/DAY 

ARGUMENT OF PERIGEE 


98.88999 

DEG 

ARGUMENT OF PERIGEE 

DERIVATIVE 

-0.28681 

DEGDAY 

MEAN ANOMALY 


20.05704 

DEG 
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1.2510762 ER 

0.79795717 
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DESCRIPTIVE SPACE ELEMENTS 



EPOCH 
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SSSS 
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71 220 000. 

0 

UT2W 

J.D.S. 
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DO 

UT2W 
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118.296797 
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PERIGEE HEIGHT 670.1946 KM 

416.4396 


ST MI 

APOGEE HEIGHT 2532.6167 KM 

1573.6951 


ST MI 

NQOE— SUN- ANGLE , R AA NODE MINUS RA SUN 

15.37797 


DEG 

LONGITUDE OF ASCENDING NODE 

161.61571 


DEG 

LATITUOE OF P ER I GEE , GEOCE NTRI C 

77.138428 


DEG 

VELOCITY AT PERIGEE 7.946854 KM/SC 
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ST MI/HR 

VELOCITY AT APOGEE 6.28590C KM/SC 
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PREDICTION 

SPACE ELEMENTS 


EPOCH 


T- 1 

T- 2 

T- 3 



YRMODYHRMMSS. SSSS 

YRMODYHRMMSS. SSSS 

YRMODYHRMMSS. SSSS 

CAL DAT UT2W 


719223 0. 0 

710302 0. 0 

710359 C. 0 

J.D.S. UT2W 


4906 .OOOCOOOOO 

4913.000000000 

4920.000000000 

PERIOD 

MIN 

118.288893 

118.288484 

118.288042 

PERIOD DER MD/D 

-". 0405 3 

-0.04385 

-0.04674 

ECCENTRICITY 


0.11671034 

0.11668890 

0.11660359 

..I NCI INATION 

DEG 

80.667821 

80.667843 

80.667935 

RA ASC NODE 

DEG 

345.392596 

340.103199 

334.813777 

ARG PERIGEE 

DEG 

92.867035 

78.813560 

64.752881 

MEAN ANOM 

DEG 

207.520943 
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2.610200 
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ER 
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1.2510722 

1.2510694 
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OSCULATING SPACE ELEMENTS 


EPOCH 



YRMODYHRMMSS.SSSS 


CAL OAT 


71 220 COO. 0 

UT2W 

J • D . S 

• 


4903. 00 

UT2W 

OSCULATING SPACE ELEMENTS 



PERIOD 


118.116753 

Ml N 

ECCENTRICITY 
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INCL IN AT ION 


80.66564 

DEG 
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ASCENSION OF 

ASCENDING NODE 

347.65290 

DEG 

ARGUMENT OF PERIGEE 

98.50309 

DEG 

MEAN 

ANOMALY 


20.39206 

DEG 

OSCULATING CARTESIAN UUANTIES 



X 

-0.37110174 
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ER 

Y 
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ER 

Z 

0.58105523 

DECAMEGAME TERS 
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ER 

vx 

-0.57792410 
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ER/CUT 

VY 
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ER/CUT 

VZ 

-0.35186586 

DMM/CENTIDAY 
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ER/CUT 

GM 
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MEAN ELEMENTS OF MODIFIED BROUWER 

TYPE 


EPOCH 

YRMOD YHRMMSS . 
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CAL DAT 

71 220 C 00. 

* 

UT2W 

J « 0 • s • 
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00 
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ER 
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EPOCH YRMODYHRMMSS.SSSS 
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o 

o 

0 UT2W 
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00 
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04 
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04 
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04 
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03 
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05 
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04 
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04 
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05 
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04 
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05 
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05 
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05 
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05 


F igure 4e 


191 
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Figure 4d 


E QUA T 0^ CROSSINGS 
NORTHBOUND 
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11306 

2 2.44 
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11307 
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11308 
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11309 
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11310 
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11311 
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11312 
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17.70 

11313 
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47.43 

11314 
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11315 
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11316 
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11317 

2344.21 
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11319 
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11321 
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11323 
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11326 
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11327 
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11328 
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11330 

122.64 

1 92 .82 

11331 

320.98 

222.55 
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11333 
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11334 
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11335 
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341.46 

11336 
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11.19 

11337 

1511.03 

40.92 

11338 
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11339 
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11341 
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11343 
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2244.48 

156.57 

28 FEB 71 








11354 

C42.32 

1 86.30 

11355 

241.16 

216.02 

11356 

439.50 

245.75 

11357 

637. 84 

275.48 

113 58 

836.1 8 

305.21 

11359 

1034.52 
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16 9.61 

61 .86 

11375 

18 7.95 
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0 2.97 

179.77 

11379 

2 1.31 

209.50 

11383 
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ONE ORBIT EPHEMERIS 
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820.4* 

SN 

60 
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92.13 

49.80 

2466.4* 

SN-6D 

96 .49 

40 ,67 

2362.3* 

SN-5 0 

100.57 
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1922.5* 

SN-20 

111.70 

31 .47 
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8ROUWER CAP PARAMETERS 
NO OAF PARAMETERS COMPUTED IN THIS RUN 
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